Methods for preparing entangled quantum states

ABSTRACT

Various embodiments of the present invention are directed to methods for generating an entangled state of qubits. In one embodiment of the present invention, a method for preparing an entangled state of qubits comprises providing a probe and N non-interacting qubits, each qubit comprises a linear superposition of two basis states. The probe is transmitted into an interaction region that separately couples the probe to each of the qubits and produces a number of different probes. A linear superposition of states is output from the interaction region, each state in the linear superposition of states comprises a tensor product of entangled basis states and one of the different probes. The linear superposition of states is projected into one of the entangled states by measuring the state of the probe.

TECHNICAL FIELD

The present invention relates to quantum computation and quantum information, and, in particular, to methods for preparing various entangled states of qubits.

BACKGROUND OF THE INVENTION

A bit is a basic unit of information used by conventional computational systems to process information and store information in information-storage devices, such as magnetic and optical discs. The bit is equivalent to a choice between two mutually exclusive alternatives, such as “on” and “off,” and is typically represented by the numbers 0 or 1. Information encoded in bits is called “classical information.” In recent years, advancements in the field of physics have given rise to methods of encoding information in discrete, or continuous, states of quantum systems, including electrons, atoms, and photons of electromagnetic radiation. Information encoded in the discrete states of a quantum system is called “quantum information.” An elementary quantum system has two discrete states and is called a “qubit.” The qubit “basis states” are represented by

and

and are used to represent the bits 0 and 1, respectively. However, unlike the systems used to realize bits in classical information, such a quantum system can be in the state

the state

or in a state that simultaneously comprises both

and

These qubit states are represented by a linear superposition of states:

The parameters α and β are complex-valued coefficients satisfying the condition: |α|²+|β|²=1 where |α|² is the probability of measuring the state

and |β|² is the probability of measuring the state

A qubit can exist in any one of an infinite number of linear superpositions until the qubit is measured. When the qubit is measured in the computational basis

and

the qubit is projected into either the state

or the state

The infinite number of qubit-linear superpositions can be geometrically represented by a unit-radius, three-dimensional sphere called a “Bloch sphere”: $\left. \psi \right\rangle = {{{\cos\left( \frac{\theta}{2} \right)}\left. 0 \right\rangle} + {{\mathbb{e}}^{{\mathbb{i}}\quad\phi}\quad{\sin\left( \frac{\theta}{2} \right)}\left. 1 \right\rangle}}$ where −π/2<θ<π/2 and 0<φ≦π. FIGS. 1A-1C illustrate a Bloch sphere representation of a qubit. In FIG. 1A, lines 101-103 are the orthogonal x, y, and z Cartesian coordinate axes, respectively, and the Bloch sphere 106 is centered at the origin. There are an infinite number of points on the Bloch sphere 106, each point representing a unique linear superposition of the qubit

For example, a point 108 on the Bloch sphere 106 represents a qubit

comprising in part the state

and in part the state

However, once the state of the qubit

is measured in the computational basis

and

the qubit

is projected into the state

110, in FIG. 1B, or the state

112, in FIG. 1C.

Two or more quantum systems can be used to encode bit strings. For example, the four, two-bits strings “00,” “01,” “10,” and “11” can be correspondingly encoded in the two-qubit product states

and

where the subscript “1” represents a first qubit system, and the subscript “2” represents a second qubit system. However, the first qubit system and the second qubit system can exist simultaneously in two basis states that are represented by a linear superposition of the product states as follows: $\begin{matrix} {\left. \psi \right\rangle = {\frac{1}{2^{1/2}}\left( {\left. 0 \right\rangle_{1} + \left. 1 \right\rangle_{1}} \right)\frac{1}{2^{1/2}}\left( {\left. 0 \right\rangle_{2} + \left. 1 \right\rangle_{2}} \right)}} \\ {= {\frac{1}{2}\left\lbrack {{\left. 0 \right\rangle_{1}\left. 0 \right\rangle_{2}} + {\left. 0 \right\rangle_{1}\left. 1 \right\rangle_{2}} + {\left. 1 \right\rangle_{1}\left. 0 \right\rangle_{2}} + {\left. 1 \right\rangle_{1}\left. 1 \right\rangle_{2}}} \right\rbrack}} \end{matrix}$ The state

indicates that by squaring the coefficient ½ there is a ¼ probability of measuring each of the product states

and

when the two qubits are measured separately, each in their computation basis. Certain linear superpositions of the product states, called “entangled states,” can be used in quantum computing and to process and transmit quantum-information. Quantum entanglement is a quantum mechanical property in which the states of two or more quantum systems are linked to one another, even though the quantum systems may be spatially separated. Such entangled states cannot be written as a simple product of a state for each system. The following linear superpositions, called “the Bell states”: ${\left. \psi^{-} \right\rangle_{12} = {\frac{1}{\sqrt{2}}\left( {{\left. 0 \right\rangle_{1}\left. 1 \right\rangle_{2}} - {\left. 1 \right\rangle_{1}\left. 0 \right\rangle_{2}}} \right)}},{\left. \psi^{+} \right\rangle_{12} = {\frac{1}{\sqrt{2}}\left( {{\left. 0 \right\rangle_{1}\left. 1 \right\rangle_{2}} + {\left. 1 \right\rangle_{1}\left. 0 \right\rangle_{2}}} \right)}},{\left. \phi^{-} \right\rangle_{12} = {\frac{1}{\sqrt{2}}\left( {{\left. 0 \right\rangle_{1}\left. 0 \right\rangle_{2}} - {\left. 1 \right\rangle_{1}\left. 1 \right\rangle_{2}}} \right)}},{and}$ $\left. \phi^{+} \right\rangle_{12} = {\frac{1}{\sqrt{2}}\left( {{\left. 0 \right\rangle_{1}\left. 0 \right\rangle_{2}} + {\left. 1 \right\rangle_{1}\left. 1 \right\rangle_{2}}} \right)}$ are examples of entangled states. Consider a first qubit system and a second qubit system that have both been prepared in the Bell state

The square of the coefficient 1/√{square root over (2)} indicates that when a measurement is performed to determine the state of the first and second qubit systems, there is a ½ probability of obtaining the result

and a ½ probability of obtaining the result

Suppose that after the quantum systems have been spatially separated, the first qubit system is measured and determined to be in the state

Quantum entanglement ensures that the second qubit system is measured in the state

in spite of the fact that the two qubit systems are spatially separated and measured at different times.

Entangled qubit states have a number of different and useful quantum-enhanced applications, such as quantum metrology, quantum cryptography, quantum communication, and quantum teleportation. For the sake of simplicity, quantum teleportation is described below as an example of a quantum-enhanced application. Quantum teleportation can be used to transmit quantum information in the absence of a quantum communications channel linking the sender of the quantum information to the recipient of the quantum information. FIG. 2 illustrates an example of quantum teleportation. In FIG. 2, quantum communications channels are represented by directional arrows, such as directional arrow 202, and a classical communications channel is represented by a dashed-line directional arrow 204. Bob receives a qubit

that he needs to transmit to Alice, but Bob does not know the value of the parameters α and β, and Bob can only transmit classical information over the communications channel 204. Bob can perform a measurement on the qubit, but by the properties of quantum mechanics, a measurement irrevocably destroys most of the qubit information rendering any information that Alice may choose to obtain from

irretrievable. In order for Alice to receive all of the information contained in the qubit

both Bob and Alice agree in advance to share a Bell state generated by an entangled state source 206:

=1/√{square root over (2)}(|0>_(A)

_(B)) where the subscript “A” identifies qubit basis states transmitted to Alice, and the subscript “B” identifies qubit basis states transmitted to Bob.

The overall state the system in FIG. 2 is:

=1/√{square root over (2)}(α|0>_(C)+β|1>_(C))(|0>_(A)|1>_(B)−

0>_(B)) where the qubit

is represented by

The state

can be rewritten in terms of the four Bell states as follows: $\left. \Theta \right\rangle_{ABC} = {\frac{1}{2}\left\lbrack {{\left. \psi^{-} \right\rangle_{BC}\left( {{\alpha\left. 0 \right\rangle_{A}} - {\beta\left. 1 \right\rangle_{A}}} \right)} + {\left. \psi^{+} \right\rangle_{BC}\left( {{{- \alpha}\left. 0 \right\rangle_{A}} - {\beta\left. 1 \right\rangle_{A}}} \right)} + {\left. \phi^{-} \right\rangle_{BC}\left( {{\beta\left. 0 \right\rangle_{A}} - {\alpha\left. 1 \right\rangle_{A}}} \right)} + {\left. \phi^{+} \right\rangle_{BC}\left( {{{- \beta}\left. 0 \right\rangle_{A}} - {\alpha\left. 1 \right\rangle_{A}}} \right)}} \right\rbrack}$ The state

reveals that the Bell states are entangled with the qubits identified by the subscript “A.” Bob is in possession of the Bell states identified by the subscript “BC,” and Alice is in possession of the qubits identified by the subscript “A,” but Alice does not know which of the four qubit states she possesses. Bob and Alice both agree in advance that the strings “00,” “01,” “10,” and “11” correspond to the entangled states

and

so that when Bob performs a Bell state measurement to determine the Bell states in his possession, he can immediately transmit to Alice the corresponding two-bit string over the communications channel 204. As a result, Alice knows which qubit state she possesses. For example, suppose Bob performs a Bells state measurement that outputs the state

Bob's measurement projects the state

into the state

(α|0>_(A)−β

). Bob then transmits the string “01” over the communications channel 204 to Alice. Quantum entanglement ensures that Alice knows with certainty that she possesses the qubit state

which is equivalent to the original qubit state as the overall phase is unimportant. For the other Bell state measurement outcomes, all of which occur with probability ¼, Alice performs operations on the qubit in order to transform the state into the original unknown state supplied by Bob.

Methods used to generate entangled qubits often employ direct interactions between the qubits being entangled, are limited to either matter-based qubits or photon-based qubits, or necessitate performing numerous measurements. As a result, physicists, computer scientists, and users of quantum information have recognized a need for new methods that can be used to generate entangled qubits in both matter-based and photon-based qubits using a single measurement and indirect interactions between qubits.

SUMMARY OF THE INVENTION

Various embodiments of the present invention are directed to methods for generating an entangled state of qubits. In one embodiment of the present invention, a method for preparing an entangled state of qubits comprises providing a probe and N non-interacting qubits, each qubit comprises a linear superposition of two basis states. The probe is transmitted into an interaction region that separately couples the probe to each of the qubits and produces a number of different probes. A linear superposition of states is output from the interaction region, each state in the linear superposition of states comprises a tensor product of entangled basis states and one of the different probes. The linear superposition of states is projected into one of the entangled states by measuring the state of the probe.

BRIEF DESCRIPTION OF THE DRAWINGS

FIGS. 1A-1C illustrate a Bloch sphere representation of a qubit.

FIG. 2 illustrates an example of quantum teleportation.

FIG. 3 illustrates a cubic cavity.

FIG. 4 illustrates a three-dimensional coordinate system with two independent polarization vectors and a normalized wavevector as basis vectors.

FIG. 5 illustrates a representation of electric and magnetic field components of electromagnetic radiation in the coordinate system shown in FIG. 4.

FIG. 6 is an energy-level diagram of quantized electromagnetic radiation.

FIG. 7 is a phase space diagram of a number state.

FIG. 8 illustrates a probability distribution function associated with a pulse output from a source and transmitted to a detector.

FIG. 9 is a phase-space diagram of a coherent state.

FIG. 10 is a plot of a probability distribution function associated with a coherent state.

FIG. 11 illustrates a control-flow diagram that describes one of many embodiments of the present invention for generating an entangled state of non-interacting qubits via interaction with a coherent state.

FIGS. 12A-12B are phase-space diagrams that represent a phase shift and a translation of a probe, respectively.

FIG. 13 illustrates a first implementation of the method described in FIG. 11 that represents an embodiment of the present invention.

FIG. 14A is a phase-space diagram of probes entangled with basis states output from a probe-basis-state interaction region of the implementation in FIG. 13 that represents an embodiment of the present invention.

FIG. 14B is a plot of three probability distributions associated with the phase-space diagram in FIG. 13A that represents one of many embodiments of the present invention.

FIGS. 15A-15B illustrate corrections that remove a phase shift from entangled states output from the implementation shown in FIG. 13 and represents an embodiment of the present invention.

FIG. 16 illustrates a second implementation of the method described to FIG. 11 that represents an embodiment of the present invention.

FIG. 17A is a phase-space diagram of probes entangled with basis states output from a probe-basis-state interaction region of the implementation in FIG. 15 that represents an embodiment of the present invention.

FIG. 17B is a phase-space diagram of rotated probes entangled with basis states output from a probe-basis-state interaction region of the implementation in FIG. 15 that represents an embodiment of the present invention.

FIG. 17C is a plot of three probability distributions associated with the phase-space diagram of the output state in FIG. 17B that represents one of many embodiments of the present invention.

FIGS. 18A-18B illustrate corrections that remove a phase shift from entangled states output from the implementation shown in FIG. 16 and represents an embodiment of the present invention.

FIG. 19 illustrates a third implementation of the method described to FIG. 11 that represents an embodiment of the present invention.

FIG. 20A is a phase-space diagram of probes entangled with basis states output from a probe-basis-state interaction region of the implementation in FIG. 19 that represents an embodiment of the present invention.

FIG. 20B is a phase-space diagram of the probes in FIG. 20A after application of a translation operation that represents one of many embodiments of the present invention.

FIGS. 21A-21B illustrate corrections that remove a phase shift from entangled states output from the implementation shown in FIG. 19 and represents an embodiment of the present invention.

FIG. 22A is a phase-space diagram of probes entangled with basis states output from a probe-basis-state interaction region of the implementation in FIG. 19 that represents an embodiment of the present invention.

FIG. 22B is a phase-space diagram of the probes in FIG. 22A after application of a translation operation that represents one of many embodiments of the present invention.

FIG. 23 illustrates an interaction region comprising a single interaction medium that represents one of many embodiments of the present invention.

DETAILED DESCRIPTION OF THE INVENTION

Various embodiments of the present invention are directed to methods for generating entangled states of qubits. In order to assist in understanding descriptions of various embodiments of the present invention, an overview of quantum mechanics is provided below, in a first subsection. In a second subsection, an overview of electromagnetic radiation and quantum optics is provided. In a third subsection, an overview of coherent states is provided. In a fourth subsection, an overview of quantum entanglement is provided. Finally, in a fifth subsection, various method embodiments of the present invention are described.

Overview of Quantum Mechanics

Embodiments of the present invention employ concepts in quantum mechanics. The textbook “Modern Quantum Mechanics Revised Edition,” J. J. Sakurai, Addison Wesley Publishing Company, New York, 1994, is one reference for the field of quantum mechanics. In this subsection, topics in quantum mechanics that relate to embodiments of the present invention are described. Additional details can be obtained from the above-referenced textbook, or from many other textbooks, papers, and journal articles related to quantum mechanics.

Quantum mechanics models the observed behavior of systems at the atomic and subatomic levels, which comprise photons, electrons, atoms, and molecules. Quantum systems exist in discrete states that are characterized by discrete measurable quantities. A state of a quantum system is represented by a ket and is denoted

where Ψ is a label that represents a state of a quantum system. For example, an electron has two intrinsic spin-angular-momentum states that correspond to two measurable spin-angular-momentum values h/2 and −h/2, where h is approximately 1.0546×10⁻³⁴ Js. The spin state that corresponds to the spin-angular momentum h/2 is referred to as “spin up” and is denoted

and the spin state that corresponds to the spin angular momentum −h/2 is referred to as “spin down” and is denoted

Various different labels can be assigned to various different quantum states. For example, the spin up and spin down states

and

can also be represented by the kets

and

respectively. Also, a single label can be used to represent different states in entirely different quantum systems. For example, the ket

can represent a first quantized vibrational level of a diatomic molecule and can be used to represent a single photon, as described below, in a following subsection.

A measurement employed to determine a measurable quantity of a quantum system, such as the spin angular momentum of an electron, is represented by an operator

where the symbol

denotes an operator. In general, an operator operates on a ket from the left as follows:

where

is a ket representing an observed quantum state. Typically, an operator

is associated with a set of states called “eigenstates.” An eigenstate is represented as

with the following property:

where

i is a non-negative integer, and

ψ_(i) is a real value, called an “eigenvalue,” that corresponds to a discrete measurable quantity that is observed when the quantum system is in the eigenstate

For example, a measurement employed to determine the spin angular momentum of an electron is represented by Ŝ_(z), and the eigenvalue-eigenstate representations of observed spin-angular-momentum values are: ${{{\hat{S}}_{z}\left. \uparrow \right\rangle} = {\frac{\hslash}{2}\left. \uparrow \right\rangle}},{and}$ ${{\hat{S}}_{z}\left. \downarrow \right\rangle} = {{- \frac{\hslash}{2}}\left. \downarrow \right\rangle}$

The eigenstates are basis vectors of a complex vector space called a “Hilbert space,” and the number of eigenstates is the dimension of the Hilbert space. For example, a Hilbert space of an electron is two-dimensional, with eigenstates

and

A Hilbert space with N eigenstates

is N-dimensional, and any state

in the Hilbert space can be written as a linear superposition of the eigenstates as follows: $\left. \Psi \right\rangle = {\sum\limits_{i = 1}^{N}{c_{i}\left. \psi_{i} \right\rangle}}$ where c_(i) is a complex valued coefficient called the “amplitude.” A Hilbert space also includes a mathematical operation called the “inner product.” The inner product of two states

and

is represented by:

where

is called a “bra,” and represents the complex conjugate and transpose of the state

The inner product has the following property: z,57 where “*” represents the complex conjugate. The basis eigenstates of a Hilbert space are orthonormal, or in mathematical notation:

where δ_(ij) is “1” when i equals j, and 0 otherwise. For example, the inner product of the eigenstates of a single electron Hilbert space are:

and

The orthonomality property of the eigenstates of a Hilbert space can be used to determine the coefficients of the linear superposition of states

Taking the inner product of

with

gives the corresponding coefficient: $\left\langle \psi_{j} \middle| \Psi \right\rangle = {{\sum\limits_{i = 1}^{N}{c_{i}\left\langle \psi_{j} \middle| \psi_{i} \right\rangle}} = {{\sum\limits_{i = 1}^{N}{c_{i}\delta_{ij}}} = c_{j}}}$ Substituting for the coefficients in the linear superposition gives: $\left. \Psi \right\rangle = {\sum\limits_{i = 1}^{N}{\left. \psi_{i} \right\rangle\left\langle \psi_{i} \middle| \Psi \right\rangle}}$ Because

is an arbitrary ket in the Hilbert space, ${\sum\limits_{i = 1}^{N}{\left. \psi_{i} \right\rangle\left\langle \psi_{i} \right.}} = \hat{1}$ where “{circumflex over (1)}” is the identity operator. The summation is called the “completeness relation,” and the eigenstates

are said to be “complete.”

Because eigenstates of a Hilbert space are orthonormal and provide a basis for the Hilbert space, the eigenstates can be represented by orthogonal normalized column vectors and an operator can be represented by square matrix. For example, the eigenstates of a single electron Hilbert space are represented by the column vectors: ${\left. \uparrow \right\rangle\overset{.}{=}\begin{bmatrix} 1 \\ 0 \end{bmatrix}},{{{and}\quad\left. \downarrow \right\rangle}\overset{.}{=}\begin{bmatrix} 0 \\ 1 \end{bmatrix}}$ where the symbol

stands for the term “is represented by.” The complex conjugates and transposes of the eigenstates are represented by the row vectors:

and

Using the completeness relation, an operator Ô on the basis

can also be represented by: $\hat{O} = {\sum\limits_{i = 1}^{N}{\sum\limits_{j = 1}^{N}{\left. \psi_{i} \right\rangle\left\langle {\psi_{i}{\hat{O}}\psi_{j}} \right\rangle\left\langle \psi_{j} \right.}}}$ where

is a matrix element. The matrix corresponding to the operator Ô on the basis

can be represented as follows: $\hat{O}\overset{.}{=}\begin{bmatrix} \left\langle {\psi_{1}{\hat{O}}\psi_{1}} \right\rangle & \left\langle {\psi_{1}{\hat{O}}\psi_{2}} \right\rangle & \ldots & \left\langle {\psi_{1}{\hat{O}}\psi_{N}} \right\rangle \\ \left\langle {\psi_{2}{\hat{O}}\psi_{1}} \right\rangle & \left\langle {\psi_{2}{\hat{O}}\psi_{2}} \right\rangle & \quad & \vdots \\ \vdots & \quad & ⋰ & \quad \\ \left\langle {\psi_{N}{\hat{O}}\psi_{1}} \right\rangle & \ldots & \quad & \left\langle {\psi_{N}{\hat{O}}\psi_{N}} \right\rangle \end{bmatrix}$ For the operator Ô equal to the operator

the matrix representation has zero off diagonal elements, and the diagonal elements are the eigenvalues

For example, the electron spin operator can be given by: ${\hat{S}}_{z} = {\frac{\hslash}{2}{\hat{\sigma}}_{z}}$ where

The matrix representation of the electron spin operator Ŝ_(z) is given by: $\begin{matrix} {{\hat{S}}_{z}\overset{.}{=}\begin{bmatrix} \left\langle \uparrow{{\hat{S}}_{z}}\uparrow \right\rangle & \left\langle \uparrow{{\hat{S}}_{z}}\downarrow \right\rangle \\ \left\langle \downarrow{{\hat{S}}_{z}}\uparrow \right\rangle & \left\langle \downarrow{{\hat{S}}_{z}}\downarrow \right\rangle \end{bmatrix}} \\ {= {\frac{\hslash}{2}\begin{bmatrix} 1 & 0 \\ 0 & {- 1} \end{bmatrix}}} \end{matrix}$ An operator Ô that corresponds to a measurable quantity has matrix elements satisfying the condition:

=

and is said to be a “Hermitian operator.”

Prior to a measurement, a quantum system can simultaneously exist in all of the eigenstates of a corresponding Hilbert space, which is represented by the (pure state) linear superposition of states: $\left. \Psi \right\rangle = {{\sum\limits_{i = 1}^{N}{c_{i}\left. \psi_{i} \right\rangle}} = {\sum\limits_{i = 1}^{N}{\left. \psi_{i} \right\rangle\left\langle \psi_{i} \middle| \Psi \right\rangle}}}$ A measurement performed on the quantum system

projects the quantum system into one of the eigenstates. In other words, a measurement on a quantum system is essentially a filtering process that places the quantum system into one of the eigenstates in the linear superposition at the time of the measurement. For example, an electron with an unknown spin orientation prior to a measurement exists in a linear superposition of states:

A spin determination measurement Ŝ_(z) projects the linear superposition of states into either the state

or the state

at the time of the measurement.

There is a corresponding irreversible change to the state of a quantum system as a result of a measurement. Irreversibility can only be avoided when the quantum system is already in one of the quantum states before the measurement is performed. As a result, one cannot infer the prior state of a quantum system based on the outcome of a single measurement. For example, if the outcome of a spin measurement is h/2, it is not possible to determine whether the system was already in the state

or in a linear superposition of the spin states

and

at the time of the measurement.

Although it is not possible to know in advance which of the various states |ψ_(i)> a quantum system will be projected into, the probability of measuring a particular state |ψ_(i)> is given by: Probability for

=|c_(i)|²

where

is normalized, and |c_(i)|² equals c*_(i)c_(i) and gives the outcome probability. For example, prior to a spin determination measurement in the spin basis

and

consider an electron with a 1/2 probability of being in the spin state

or the spin state

The linear superposition of the electron in such as spine state prior to a spin determination measurement can be represented by: $\left. \Psi \right\rangle = {{\frac{1}{\sqrt{2}}\left. \uparrow \right\rangle} + {\frac{1}{\sqrt{2}}\left. \downarrow \right\rangle}}$

The expectation value of measurement on an ensemble of quantum systems that are described by the linear superposition of states

is mathematically represented by:

and is determined by applying the completeness relation as follows: $\begin{matrix} {\left\langle \hat{\Psi} \right\rangle = {\sum\limits_{i = 1}^{N}{\sum\limits_{j = 1}^{N}{\left\langle \Psi \middle| \psi_{i} \right\rangle\left\langle {\psi_{i}{\hat{\Psi}}\psi_{j}} \right\rangle\left\langle \psi_{j} \middle| \Psi \right\rangle}}}} \\ {= {\sum\limits_{i = 1}^{N}{\psi_{i}{\left\langle \psi_{i} \middle| \Psi \right\rangle }^{2}}}} \end{matrix}$ The expectation value represents the weighted eigenvalue average result expected from measurements on the quantum systems in the ensemble, where the initial state

of the quantum system is the same for each member of the ensemble. In other words, the linear superposition of states of each quantum system is identical prior to the measurement. In practice, such an ensemble could be realized by preparing many identical and independent quantum systems all in the same state, or by repeatedly preparing a single system in the same state. Note that the expectation value may not be the value obtained for each measurement and, therefore, is not to be confused with the eigenvalue of the measurement. For example, the expectation value of Ŝ_(z) can be any real value between the eigenvalues h/2 and −h/2, but the actual measured value of Ŝ_(z) for an electron is always either h/2 or −h/2 in each individual measurement.

A tensor product is a way of combining Hilbert spaces of different quantum systems to form Hilbert spaces that represent combined quantum systems. For example, H_(Ψ) is a Hilbert space of a first quantum system, and H_(Ξ) is a Hilbert space of a second quantum system. The Hilbert space denoted by H_(Ψ){circle around (×)}H_(Ξ) represents a combined Hilbert space, where the symbol {circle around (×)} represents a tensor product. The operators {circumflex over (Ψ)} and Ξ correspond to the Hilbert spaces H_(Ψ) and H_(Ξ), respectively, and each operates only on the corresponding eigenstates as follows:

where

represents a state in the Hilbert space H_(Ψ), and

represents a state in the Hilbert space H_(Ξ). The tensor product

can be abbreviated as

or

The spin states of two electrons in an atomic orbital is an example of a combined Hilbert space. The two electrons can either both be spin up, both be spin down, the first electron spin up and the second electron spin down, or the first electron spin down and the second electron spin up. The various tensor product representations of two spin up electrons are given by: z,99 =

=|⇑

where the subscripts 1 and 2 refer to the first and second electrons.

In quantum mechanics, there are also measurable quantities with continuous eigenvalue spectrums. The dimensionality of the corresponding Hilbert spaces are infinite and many of the properties described above for discrete quantum systems can be generalized for continuous quantum systems. A continuous eigenvalue equation is:

where ζ represents a continuous eigenvalue, and the ket

is a continuous eigenstate of the operator {circumflex over (ζ)}. For example, for an unbound particle, both position q and momentum p are continuous eigenvalues of the position and momentum operators {circumflex over (q)} and {circumflex over (p)}, respectively, and can assume any real value between −∞ and ∞.

The properties of the continuous variable ζ can be generalized as follows: ⟨ζ|ζ^(′)⟩ = δ(ζ − ζ^(′)), ∫_(−∞)^(∞)𝕕ζζ⟩⟨ζ = 1, and ⟨ζζ̂ζ^(′)⟩ = ζ^(′)δ(ζ − ζ^(′)) where ${\delta\left( {\zeta - \zeta^{\prime}} \right)} = {\lim\limits_{\Delta->0}{\frac{1}{\sqrt{2\quad\pi\quad\Delta^{2}}}{\exp\left( {- \frac{\left( {\zeta - \zeta^{\prime}} \right)^{2}}{2\quad\Delta^{2}}} \right)}}}$ A state ket for an arbitrary physical state can be expanded in terms of the states

as follows: α⟩ = ∫_(−∞)^(∞)𝕕ζζ⟩⟨ζ|α⟩ For example, consider placing in the path of a particle a detector that outputs the position of the particle when the particle is at the position q. Immediately after the measurement is taken, the system, initially in the state

is projected into the state represented by

in much the same way an arbitrary electron-spin state is projected into one of the two spin states when a spin detection measurement is performed. Other properties of the continuous variable ζ are given by: ∫_(−∞)^(∞)𝕕ζ⟨ζ|α⟩² = 1, and ⟨β|α⟩ = ∫_(−∞)^(∞)𝕕ζ⟨β|ζ⟩⟨ζ|α⟩

The momentum operator {circumflex over (p)} can also be represented by a differential operator −ih∂/∂q. As a result, both the position and momentum operators satisfy the canonical commutation relations: [q̂_(i), q̂_(j)] = 0[p̂_(i), p̂_(j)] = 0,  and[q̂_(i), q̂_(j)] = 0 = 𝕚ℏδ_(ij) where

i and j represent orthogonal coordinates, such as the Cartesian x, y, and z coordinates, and

the commutator is defined as [A, B]=AB−BA.

An Overview of Electromagnetic Radiation and Quantum Optics

In this subsection, a brief description of electromagnetic radiation and quantum optics that relates to embodiments of the present invention is described. The textbooks “Quantum Optics,” M. O. Scully and M. S. Zubairy, Cambridge University Press, Cambridge, United Kingdom, 1997, and “The Quantum Theory of Light (3^(rd) Edition),” R. Loudon, Oxford University Press, New York, 2000 are two of many references for quantum optics. Additional details can be obtained from the above-referenced textbooks, or from many other textbooks, papers, and journal articles in this field.

Quantum optics is a field of physics that relates the application of quantum mechanics to electromagnetic radiation. Electromagnetic radiation confined to a cavity with perfectly reflecting walls is quantized. Quantized electromagnetic radiation can be applied to more general unconfined optical systems, such as electromagnetic radiation propagating in free space or in an optical fiber.

Electromagnetic radiation confined to a cavity, with no free charges and currents, comprises an electric field component {right arrow over (E)}({right arrow over (r)},t) and a magnetic field component {right arrow over (B)}({right arrow over (r)},t) that are related in terms of a vector potential {right arrow over (A)}({right arrow over (r)},t) satisfying the wave equation: ${{\nabla^{2}\overset{\rightarrow}{A}} - {\frac{1}{c^{2}}\frac{\partial^{2}\overset{\rightarrow}{A}}{\partial t^{2}}}} = 0$ and the Coulomb, non-relativistic gauge condition: ∇·{right arrow over (A)}({right arrow over (r)},t)=0 where the electric and magnetic field components are determined by: ${{\overset{\rightarrow}{E}\left( {\overset{\rightarrow}{r},t} \right)} = {- \frac{\partial{\overset{\rightarrow}{A}\left( {\overset{\rightarrow}{r},t} \right)}}{\partial t}}},\quad{and}$ ${\overset{\rightarrow}{B}\left( {\overset{\rightarrow}{r},t} \right)} = {\nabla{\times {\overset{\rightarrow}{A}\left( {\overset{\rightarrow}{r},t} \right)}}}$ The electromagnetic radiation is assumed to be confined in a cubic cavity with perfectly reflecting walls, where the lengths of the walls L are much longer than the wavelengths of the electromagnetic radiation. FIG. 3 illustrates a cubic cavity 300. Orthogonal axes 302, 304, and 306 represent the x, y, and z Cartesian coordinate axes. The finite dimensional cubic cavity 300 imposes periodic boundary conditions on solutions to the wave equation. For example, in the x, y, and z-directions, plane wave solutions to the vector potential wave equation satisfy the condition: exp (i{right arrow over (k)}·{right arrow over (r)})=exp (i{right arrow over (k)}·({right arrow over (r)}+{right arrow over (L)})) where {right arrow over (L)} is (L,L,L), and

{right arrow over (k)} is called the “wavevector” with components: ${\overset{\rightarrow}{k} = {\frac{2\pi}{L}\left( {m_{x},m_{y},m_{z}} \right)}},\quad{and}$

m_(x), m_(y), and m_(z) are integers.

Each set of integers (m_(x), m_(y), m_(z)) specifies a normal mode of the electromagnetic radiation, and the magnitude of the wavevector {right arrow over (k)}, k, is equal to ω_(k)/c, where c represents the speed of light in free space and ω_(k) is the angular frequency. Note that in real life the spectrum of normal modes of an electromagnetic field is actually continuous and a discrete spectrum of normal modes suggested by the wavevector {right arrow over (k)} is an approximation to the continuous spectrum.

A vector potential solution to the wave equation above that satisfies the periodic boundary conditions is: ${A\left( {r,t} \right)} = {\sum\limits_{\overset{\rightarrow}{k},s}{{\overset{\rightarrow}{e}}_{\overset{\rightarrow}{k}s}\left( {{A_{\overset{\rightarrow}{k}s}{\mathbb{e}}^{{\mathbb{i}}{({{\overset{\rightarrow}{k} \cdot \overset{\rightarrow}{r}} - {\omega_{k}t}})}}} + {A_{\overset{\rightarrow}{k}s}^{*}{\mathbb{e}}^{- {{\mathbb{i}}{({{\overset{\rightarrow}{k} \cdot \overset{\rightarrow}{r}} - {\omega_{k}t}})}}}}} \right)}}$ where A_({right arrow over (k)}s) is a complex amplitude of the electromagnetic radiation, and {right arrow over (e_(k))}s represents two unit-length polarization vectors. The sum over {right arrow over (k)} represents the sum over the integers (m_(x), m_(y), m_(z)), and the sum over s is the sum over the two independent polarizations that are associated with each {right arrow over (k)}. The two polarization vectors are orthogonal as indicated by: {right arrow over (e_(k))}s ·{right arrow over (e_(k))}s′=δ_(ss′) and from the gauge condition given above: {right arrow over (k)}·{right arrow over (e _(k))}s=0 for both polarization directions s. The two polarization vectors {right arrow over (e_(k))}1 and {right arrow over (e_(k))}2 form a right-handed coordinate system with a normalized wavevector given by: ${{\overset{\rightarrow}{e}}_{\overset{\rightarrow}{k}1} \times {\overset{\rightarrow}{e}}_{\overset{\rightarrow}{k}2}} = {\frac{\overset{\rightarrow}{k}}{\overset{\rightarrow}{k}} = \overset{\rightarrow}{\kappa}}$ FIG. 4 illustrates a three-dimensional right-handed coordinate system with two independent polarization vectors {right arrow over (e_(k))}s and a normalized wavevector {right arrow over (k)} as basis vectors. In FIG. 4, the wavevector {right arrow over (k)} 402, and the polarization vectors, {right arrow over (e_(k))}1 404 and e_(k) 2 406, are three orthogonal unit length basis vectors of a coordinate system with coordinate axes represented by lines 408, 410, and 412, respectively.

The electric and magnetic field components of the vector potential are: ${{\overset{\rightarrow}{E}\left( {\overset{\rightarrow}{r},t} \right)} = {{\mathbb{i}}{\sum\limits_{\overset{\rightarrow}{k},s}{\omega_{k}{{\overset{\rightarrow}{e}}_{\overset{\rightarrow}{k}s}\left\lbrack {{A_{\overset{\rightarrow}{k}s}e^{{\mathbb{i}}{({{\overset{\rightarrow}{k} \cdot \overset{\rightarrow}{r}} - {\omega_{k}t}})}}} - {A_{\overset{\rightarrow}{k}s}^{*}e^{- {{\mathbb{i}}{({{\overset{\rightarrow}{k} \cdot \overset{\rightarrow}{r}} - {\omega_{k}t}})}}}}} \right\rbrack}}}}},\quad{and}$ ${\overset{\rightarrow}{B}\left( {\overset{\rightarrow}{r},t} \right)} = {\frac{\mathbb{i}}{c}{\sum\limits_{\overset{\rightarrow}{k},s}{{\omega_{k}\left( {\overset{\rightarrow}{\kappa} \times {\overset{\rightarrow}{e}}_{\overset{\rightarrow}{k}s}} \right)}\left\lbrack {{A_{\overset{\rightarrow}{k}s}e^{{\mathbb{i}}{({{\overset{\rightarrow}{k} \cdot \overset{\rightarrow}{r}} - {\omega_{k}t}})}}} - {A_{\overset{\rightarrow}{k}s}^{*}e^{- {{\mathbb{i}}{({{\overset{\rightarrow}{k} \cdot \overset{\rightarrow}{r}} - {\omega_{k}t}})}}}}} \right\rbrack}}}$ Both the electric field E( r,t) and magnetic field B( r,t) are referred to as the “classical” representation of the electric and magnetic field, are orthogonal to one another, and are both orthogonal to the wavevector k.

FIG. 5 illustrates a representation of electric and magnetic field components of electromagnetic radiation in the right-handed coordinate system shown in FIG. 4. The electromagnetic radiation is directed along the wavevector k 402 axis. The electric field component E( r,t) 502 and magnetic field component B( r,t) 504 are directed along the orthogonal polarization vectors e_(k) 1 404 and e_(k) 2 406, respectively.

The energy of the electromagnetic radiation can be determined by evaluating the Hamiltonian: $\begin{matrix} {H = {\frac{1}{2}{\int\limits_{V}{\left( {{ɛ_{0}{\overset{\rightarrow}{E} \cdot \overset{\rightarrow}{E}}} + {\frac{1}{\mu_{0}}{\overset{\rightarrow}{B} \cdot \overset{\rightarrow}{B}}}} \right){\mathbb{d}V}}}}} \\ {= {2ɛ_{0}V{\sum\limits_{\overset{\rightarrow}{k},s}{\omega_{k}^{2}A_{\overset{\rightarrow}{k}s}A_{\overset{\rightarrow}{k}s}^{*}}}}} \end{matrix}$ where ε₀ is the electric permittivity of free space,

μ₀ is the magnetic permeability of free space, and

V is the volume of the cavity.

The electric permittivity ε₀ represents the degree to which a vacuum space can store electrical potential energy under the influence of an electric field, and the magnetic permeability μ₀ represents the degree to which the vacuum modifies the flux of a magnetic field. In a non-conducting medium, the electric permittivity is further multiplied by ε, which is the degree to which the medium enhances the storage of electrical potential energy, and the magnetic permeability is further multiplied by μ, which is the degree to which the medium further enhances the flux of a magnetic field.

In order to quantize the electric field E( r,t) and magnetic field B( r,t) components, the canonical variables for position, q _(ks), and momentum, p _(ks), are introduced into the Hamiltonian by setting: $A_{\overset{\rightarrow}{k}s} = {\frac{1}{2\omega_{k}\sqrt{ɛ_{0}V}}\left( {{\omega_{k}q_{\overset{\rightarrow}{k}s}} + {{\mathbb{i}}\quad p_{\overset{\rightarrow}{k}s}}} \right)}$ As a result, the Hamiltonian for the electromagnetic radiation becomes: $H = {\frac{1}{2}{\sum\limits_{\overset{\rightarrow}{k},s}\left( {p_{\overset{\rightarrow}{k}s}^{2} + {\omega_{k}^{2}q_{\overset{\rightarrow}{k}s}^{2}}} \right)}}$ Each term in the Hamiltonian is the energy of a harmonic oscillator with vibrational mode ks, where the term p _(ks) ²/2 is the kinetic energy, and the term ω_(k)q _(ks) ²/2 is the potential energy of a harmonic oscillator with a unit mass. The Hamiltonian is quantized by replacing the position and momentum variables q _(ks) and p _(ks) with quantum mechanical position and momentum operators {circumflex over (q)} _(ks) and {circumflex over (p)} _(ks), respectively, to give the quantum Hamiltonian operator: $\hat{H} = {\frac{1}{2}{\sum\limits_{\overset{\_}{k},s}\left( {{\hat{p}}_{\overset{\_}{k}\quad s}^{2} + {\omega_{k}^{2}{\hat{q}}_{\overset{\_}{k}\quad s}^{2}}} \right)}}$ Annihilation and creation operators are defined by: ${{\hat{a}}_{\overset{\_}{k}\quad s} = {\frac{1}{\sqrt{2{\hslash\omega}_{k}}}\left( {{\omega_{k}{\hat{q}}_{\overset{\_}{k}\quad s}} + {{\mathbb{i}}\quad{\hat{p}}_{\overset{\_}{k}\quad s}}} \right)}},{and}$ ${{\hat{a}}_{\overset{\_}{k}\quad s}^{\dagger} = {\frac{1}{\sqrt{2{\hslash\omega}_{k}}}\left( {{\omega_{k}{\hat{q}}_{\overset{\_}{k}\quad s}} - {{\mathbb{i}}\quad{\hat{p}}_{\overset{\_}{k}\quad s}}} \right)}},$ and substituting the annihilation and creation operators into the quantum Hamiltonian operator gives: $\hat{H} = {\sum\limits_{\overset{\_}{k},s}{{\hslash\omega}_{k}\left( {{{\hat{a}}_{\overset{\_}{k}\quad s}^{\dagger}{\hat{a}}_{\overset{\_}{k}\quad s}} + \frac{1}{2}} \right)}}$ where â _(ks) ^(†)â _(ks) is called the “number operator” and is also denoted by {circumflex over (n)} _(ks). Using the canonical commutation relations for the position and momentum operators, the annihilation and creation operators satisfy the commutation relations given by: $\left\lbrack {{\hat{a}}_{\overset{\_}{k}\quad s},{\hat{a}}_{{\overset{\_}{k}}^{\prime}s^{\prime}}} \right\rbrack = {0 = {{\left\lbrack {{\hat{a}}_{\overset{\_}{k}\quad s}^{\dagger},{\hat{a}}_{{\overset{\_}{k}}^{\prime}s^{\prime}}^{\dagger}} \right\rbrack\left\lbrack {{\hat{a}}_{\overset{\_}{k}\quad s}^{\dagger},{\hat{a}}_{{\overset{\_}{k}}^{\prime}s^{\prime}}^{\dagger}} \right\rbrack} = {\delta_{{\overset{\_}{kk}}^{\prime}}\delta_{{ss}^{\prime}}}}}$

When the electromagnetic radiation is quantized, the amplitudes A _(ks) become operators: ${\hat{A}}_{\overset{\_}{k}\quad s} = {\sqrt{\frac{\hslash}{2\omega_{k}ɛ_{0}V}}{\hat{a}}_{\overset{\_}{k}\quad s}}$ which can be substituted into the classical electric and magnetic field equations above to obtain electric and magnetic field operators: ${{\hat{E}\left( {\overset{\_}{r},t} \right)} = {{\mathbb{i}}{\sum\limits_{\overset{\_}{k},s}{\sqrt{\frac{\hslash\omega}{2ɛ_{0}V}}{e_{\overset{\_}{k}\quad s}\left\lbrack {{{\hat{a}}_{\overset{\_}{k}\quad s}{\mathbb{e}}^{{\mathbb{i}}{({{k \cdot r} - {\omega_{k}t}})}}} - {{\hat{a}}_{\overset{\_}{k}\quad s}^{\dagger}{\mathbb{e}}^{- {{\mathbb{i}}{({{k \cdot r} - {\omega_{k}t}})}}}}} \right\rbrack}}}}},{and}$ ${\hat{B}\left( {\overset{\_}{r},t} \right)} = {\frac{\mathbb{i}}{c}{\sum\limits_{\overset{\_}{k},s}{\left( {\overset{\_}{\kappa} \times e_{\overset{\_}{k}\quad s}} \right)\sqrt{\frac{\hslash\omega}{2ɛ_{0}V}}{e_{\overset{\_}{k}\quad s}\left\lbrack {{{\hat{a}}_{\overset{\_}{k}\quad s}{\mathbb{e}}^{{\mathbb{i}}{({{k \cdot r} - {\omega_{k}t}})}}} - {{\hat{a}}_{\overset{\_}{k}\quad s}^{\dagger}{\mathbb{e}}^{- {{\mathbb{i}}{({{k \cdot r} - {\omega_{k}t}})}}}}} \right\rbrack}}}}$ Both the electric and magnetic field operators are Hermitian and represent measurable electric and magnetic fields.

Most electromagnetic radiation interactions with matter result from the electric field component rather than the magnetic field component, because the magnetic field is smaller than the electric field by the factor 1/c. As a result, the electric field alone is generally used to characterize the behavior of electromagnetic radiation and any interactions with matter, and the magnetic field component can be ignored.

Quantum computation and quantum information processing systems can be operated using a single-mode ks of electromagnetic radiation. As a result, the Hamiltonian operator for a single-mode of electromagnetic radiation reduces to: $\hat{H} = {{\hslash\omega}\left( {{{\hat{a}}^{\dagger}\hat{a}} + \frac{1}{2}} \right)}$ where â and â^(†) replace the operators â _(k) _(j) _(s) _(j) and at â _(k) _(j) _(s) _(j) ^(†) in the Hamiltonian above. The eigenstates and the corresponding energy eigenvalues of the single-mode Hamiltonian are: $\left. {\left. {\left. \hat{H} \middle| n \right\rangle = \left. {{\hslash\omega}\left( {{{\hat{a}}^{\dagger}\hat{a}} + \frac{1}{2}} \right)} \middle| n \right.} \right\rangle = \left. E_{n} \middle| n \right.} \right\rangle$ where |n> is called a “number state,” n is a nonnegative integer called a “photon number,” and E_(n) is an energy eigenvalue.

The annihilation and creation operators operate on a number state as follows: $\left. {\left. {\left. {\left. {\left. {\left. \hat{a} \middle| n \right\rangle = \left. \sqrt{n} \middle| {n - 1} \right.} \right\rangle,\left. {\hat{a}}^{\dagger} \middle| n \right.} \right\rangle = \left. \sqrt{n + 1} \middle| {n + 1} \right.} \right\rangle,\left. {{and}\hat{n}} \middle| n \right.} \right\rangle = \left. n \middle| n \right.} \right\rangle$ where {circumflex over (n)} represents the operator â^(†)â and is called the “number operator.” The number states can be generated by repeated application of the annihilation and creation operators to the number states. For example, repeated application of the annihilation operator to a number state lowers the photon number: $\left. {{{\left. 0 \right\rangle\frac{{\hat{a}}^{n}}{\sqrt{n!}}}}n} \right\rangle$ where |0> is called the “vacuum state,” which represents the lowest energy state of the electromagnetic radiation. Beginning with the vacuum state, and repeatedly applying the creation operator gives: $\left. n \right\rangle = {\frac{\left( {\hat{a}}^{\dagger} \right)^{n}}{\sqrt{n!}}\left. 0 \right\rangle}$ The number states are orthogonal and form a compete set represented by: ${\left\langle {n^{\prime}❘n} \right\rangle = \delta_{n^{\prime}n}},{{{and}\quad{\sum\limits_{n = 0}^{\infty}\quad{\left. n \right\rangle\left\langle n \right.}}} = 1}$

In general, the energy eigenvalue equation associated with a number state |n> is: ${\hat{H}\left. n \right\rangle} = {{{{\hslash\omega}\left( {n + \frac{1}{2}} \right)}\left. n \right\rangle} = {E_{n}\left. n \right\rangle}}$ Applying the annihilation and creation operators to the energy eigenvalue equation gives: ${{\hat{H}\left( {\hat{a}\left. n \right\rangle} \right)} = {{{{\hslash\omega}\left( {n - \frac{1}{2}} \right)}\left. {n - 1} \right\rangle} = {\left( {E_{n} - {\hslash\omega}} \right)\left. {n - 1} \right\rangle}}},{{{and}\quad{\hat{H}\left( {{\hat{a}}^{\dagger}\left. n \right\rangle} \right)}} = {{{{\hslash\omega}\left( {n + \frac{3}{2}} \right)}\left. {n + 1} \right\rangle} = {\left( {E_{n} + {\hslash\omega}} \right)\left. {n + 1} \right\rangle}}}$ which shows that the energy levels of electromagnetic radiation are equally spaced by a quantum of energy hω. In other words, the excitations of electromagnetic radiation occur in discrete amounts of energy hω called “photons.” The photon number n refers to the number of photons hω comprising the electromagnetic radiation.

FIG. 6 is an energy level diagram of quantized electromagnetic radiation. Horizontal lines, such as horizontal line 602, represent energy levels of electromagnetic radiation. Energy level 604 is the lowest energy level, which corresponds to the vacuum state

The energy of the vacuum state is hω/2 or 1/2 the energy of a single photon. Higher energy levels of electromagnetic radiation are each separated by the same quantum of energy hω. For example, the energy level 606 represents electromagnetic radiation with a total electromagnetic energy of 5 hω/2, which can be thought of as the energy of two photons plus the vacuum state energy hω/2. The annihilation operator corresponds to removal of a photon from the electromagnetic radiation, and the creation operator corresponds to addition of a photon to the electromagnetic radiation. For example, the annihilation operator â represents an electromagnetic-radiation transition 608 from the state

602 to the lower energy state

610. The transition 608 is achieved by giving up a photon to the surroundings. By contrast, the creation operator â^(†) represents an electromagnetic-radiation transition 612 from the state

602 to the higher energy state

614. The transition 612 is achieved by accepting a photon from the surroundings. Note that typically the surroundings can be an atom, a quantum dot, or any other system that couples to the field through a dipole interaction. Loss or absorption of a photon will involve a simultaneous excitation of the surrounding system and creation or emission of a photon will involve a corresponding de-excitation of the surrounding system.

Both the creation and annihilation operators are not Hermitian. As a result, the operators â and â^(†) cannot represent measurable quantities. However, the annihilation and creation operators can be used to construct the following Hermitian quadrature operators: ${\hat{X} = {{\sqrt{\frac{m\omega}{2\hslash}\quad}\hat{q}} = {\frac{1}{2}\left( {\hat{a} + {\hat{a}}^{\dagger}} \right)}}},{{{and}\quad\hat{Y}} = {\frac{\hat{p}}{\sqrt{2{m\hslash\omega}}} = {\frac{1}{2}{i\left( {{\hat{a}}^{\dagger} - \hat{a}} \right)}}}}$ The quadrature operators are essentially dimensionless position and momentum operators and are associated with the electric field amplitudes oscillating out of phase with each other by 90°. The energy eigenvalue can be rewritten in terms of the quadrature operators as: Ĥ=hω({circumflex over (X)} ² +Ŷ ²)

=hω(n+½)

The number states have the quadrature-operator eigenvalue property: ({circumflex over (X)} ² +Ŷ ²)

=(n+½)

and the number states have identical properties for the {circumflex over (X)} and Ŷ quadrature operators. For example, the quadrature-operator-expectation values:

{circumflex over (X)}=<n|Ŷ=0 The quadrature operators can be used to construct a phase-space diagram of the number states.

FIG. 7 is a phase space diagram of a number state. Horizontal line 702 is the X quadrature axis, and vertical line 704 is the Y quadrature axis. The origin 706 corresponds to the vacuum state

and circle 708 is a qualitative representation of a number state

where the radius of the circle 708 is the number of photons n. The number state

can be located anywhere on the circle 708, which indicates that uncertainty in the number of photons is zero, but the phase angle associated with the number state is entirely undefined, so measurement will give all possible outcomes with equal probability.

The number states also have the property: Σ=

E (t)|n>=0 where Σ is called a coherent signal. The zero valued coherent signal of a photon state is consistent with the sinusoidal variation of the electric field with time at a fixed observation point.

Photons can be generated by a photon source and transmitted through free space or in an optical fiber. The photon source can be a pulsed laser that generates a single pulse or a train of pulses, each pulse containing one or more photons that all have the same optical properties, such as wavelength and direction. Photons with the same optical properties are called “coherent.” However, the source, the detector, and a medium, such as an optical fiber, separating the source from the detector do not define an optical cavity. The source and the detector are parts of a continuous unidirectional flow of optical energy with no significant reflection or recycling of the optical energy. A pulse transmitted through free space or an optical fiber is described by a wavepacket that can be represented by a time-dependent, Gaussian-shaped function given by: ${\xi\quad(t)} = {\left( \frac{2\Delta^{2}}{\pi} \right)^{\frac{1}{4}}\exp\left\{ {{{- {\mathbb{i}\omega}_{0}}t} - {\Delta^{2}\left( {t_{0} - t} \right)}^{2}} \right\}}$ where ω₀ is the central frequency of the pulse spectrum, t is time, t₀ is the time at which the peak of the wavepacket is located at a distance z₀ from the photon source, and Δ² is the variance of the intensity spectrum. The time t₀ can be determined by z₀/v, where v is the velocity of the pulse traveling through free space or in an optical fiber.

The wavepacket ξ(t) is the amplitude of the pulse, and |ξ(t)|² is a photodetection probability density function of the pulse, where the photodetection probability density function |ξ(t)² satisfies the normalization condition: ∫_(−∞)^(∞)  𝕕tξ  (t)² = 1 The probability of photodetection of a photon in the time interval (t₁,t₂) at a distance z₀ from the photon source is given by: Probability  of  (t₁ < t₂) = ∫_(t₁)^(t₂)  𝕕tξ  (t)²

FIG. 8 illustrates a probability distribution associated with a pulse output from a source 802 and transmitted in an optical fiber 804 to a detector 806. The horizontal line 808 represents the distance zo the photon travels from the source 802 to the detector 806, and the horizontal line 810 is a time axis. The curve 812 represents the photodetection probability density function |ξ(t)². In FIG. 8, the photodetection probability density function |ξ(t)² 812 is centered at the time t₀, which corresponds to the time a pulse takes to travel the distance z₀. An area under the curve 812 represents the probability of detecting the pulse within a particular time period. For example, hash-marked region 814 represents the probability of detecting the photon within the time period t₁<t₀<t₂. Time period 816 is called a “time bin” and corresponds to the time period within which the photon is detected at the detector 806.

The time dependent creation operators can be used to generate a photon wavepacket creation operator as follows: â_(ξ)^(†) = ∫_(−∞)^(∞)  𝕕t  ξ  (t)â^(†)(t) The creation operator can be used to construct continuous-mode number states that represent photons transmitted through free space or in an optical fiber as follows: $\left. n_{\xi} \right\rangle = {\frac{\left( {\hat{a}}_{\xi}^{\dagger} \right)^{n}}{\sqrt{n!}}\left. 0 \right\rangle}$ where |0> is the continuous-mode vacuum state. The continuous-mode number states satisfy the following same conditions: n̂n_(ξ)⟩ = nn_(ξ)⟩ ⟨n_(ξ)^(′)|n_(ξ)⟩ = δ_(n^(′)n), and ${\sum\limits_{n_{\xi} = 0}^{\infty}{\left. n_{\xi} \right\rangle\left\langle n_{\xi} \right.}} = 1$ As a result, the subscript ξ used to identify continuous-mode number states can be dropped. Note that the wavepacket constructed photon is not an eigenstate of any Hamiltonian.

An Overview of Coherent States

The most common kind of single-mode states are linear superpositions of the number states. There are a number of different possible linear superpositions of the number states, but the coherent state: $\left. \alpha \right\rangle = {{\exp\left( {{- \frac{1}{2}}{\alpha }^{2}} \right)}{\sum\limits_{n = 0}^{\infty}{\frac{\alpha^{n}}{\sqrt{n!}}\left. n \right\rangle}}}$ is a linear superposition of the number states used in many applications of quantized electromagnetic radiation. The coherent states are eigenstates of the annihilation operator:

where taking the complex conjugate gives:

However, the coherent state

is not an eigenstate of the creation operator at because the summation over a cannot be rearranged to give the coherent state from

The coherent state expectation value for the number operator:

{circumflex over (n)}=|α| ² indicates that |α|² is the mean number of photons. The probability of detecting n photons in a measurement of the number of photons is a Poisson distribution: $P_{n} = {{\left\langle n \middle| \alpha \right\rangle }^{2} = {{\exp\left( {- {\alpha }^{2}} \right)}\frac{{\alpha }^{2}}{n!}}}$ The Poisson distribution approaches a Gaussian distribution for large values of |α|².

The coherent state is a quantum state whose properties most closely resemble a classical electromagnetic wave of stable amplitude and fixed phase. For example, the electric field operator corresponding to an electric field propagating in the z direction, with the mode subscripts k and s removed, is: $\begin{matrix} {{\hat{E}(\Omega)} = {\frac{1}{2}\left( {{\hat{a}{\mathbb{e}}^{{- {\mathbb{i}}}\quad\Omega}} + {{\hat{a}}^{\dagger}{\mathbb{e}}^{{\mathbb{i}}\quad\Omega}}} \right)}} \\ {= {{\hat{X}\cos\quad\Omega} + {\hat{Y}\sin\quad\Omega}}} \end{matrix}$ where the time t and displacement z are contained in the phase angle: ${\Omega\left( {z,t} \right)} = {{\omega\quad t} - {kz} - \frac{\pi}{2}}$ and the electric field is measured in units of √{square root over (hω/2ε₀V)}.

The coherent state is a nearly classical state because it gives the correct sinusoidal form for the electric field expectation value or coherent signal: $\begin{matrix} {\Sigma = \left\langle {\alpha{{\hat{E}(\Omega)}}\alpha} \right\rangle} \\ {= {{\alpha }{\cos\left( {\Omega - \phi} \right)}}} \end{matrix}$ where α=|α|e^(iφ), and

φ is the mean phase angle of the coherent state excitation of the mode.

FIG. 9 is a phase-space diagram of a coherent state. Horizontal axis 902 is the X quadrature axis, and vertical axis 904 is the Y quadrature axis. Disk 906 represents an area of uncertainty associated with the phase and average number of photons of a coherent state

the fluctuations being equal in all directions of phase space. Directional arrow 908 represents the amplitude |α|, or the mean number of photons, and (Ω-φ) is the angle above the X quadrature axis. The uncertainty disk is resolved into a photon-number uncertainty and phase uncertainty, as indicated by directional arrows 910 and 912, respectively. Directional arrow 910 represents uncertainty in the phase angle φ, which diminishes as |α| increases according to: ${\Delta\quad\phi} = \frac{1}{2\quad{\alpha }}$ Because the phase uncertainty is inversely proportional to |α|, the coherent state becomes better defined as the average number of photons is increased. Directional arrow 912 represents the photon-number uncertainty: Δn=|α| The diameter of the uncertainty disk is: ${\Delta\quad\phi\quad\Delta\quad n} = \frac{1}{2}$ A measurement on the coherent state outputs the coherent signal Σ, which is represented by projecting the center of the uncertainty disk 906 onto the X-quadrature axis 902.

The angle Ω is a property of the measurement that can be set equal to zero by the experimentalist, which gives an X-quadrature-expectation value:

=|α|cos φ=K A homodyne detection measurement outputs the X-quadrature-expectation value K as a function of the phase angle φ and the amplitude |α|. The quantity measured is the difference between the numbers of photons arriving at two different photodetectors during a period of time. The homodyne detector measures photon numbers, or photon counts, and the effect is to produce measurements proportional to an electric field quadrature, enabled through the measurement of detector currents.

A probability amplitude associated with homodyne detection of the coherent state is given by a Gaussian function: ${f\left( {x,\beta} \right)} = {\frac{1}{\left( {2\quad\pi} \right)^{1/4}}{\exp\left( {- \frac{\left( {x - {2\quad\beta}} \right)^{2}}{4}} \right)}}$ where β is equal to α cos φ, and

x is the value signal output from the homodyne detection.

Squaring the amplitude f (x,β) gives the probability distribution of a homodyne measurement. FIG. 10 is a plot of a probability distribution 1002 associated with the coherent state representation shown in FIG. 9. Horizontal axis 1004 corresponds to X quadrature values, and vertical axis 1006 is a probability axis. The probability distribution 1002 is centered about 2α cos φ. A coherent signal output by a homodyne measurement falls within the region 1008. The probability that the coherent signal falls within a particular interval is associated with an area under the probability distribution 1002. For example, hash-marked region 1010 represents the probability that the homodyne detection gives a value in the interval 1012.

An Overview of Quantum Entanglement

A quantum system comprising a first quantum subsystem and a second quantum subsystem has a Hilbert space H_(A){circle around (×)}H_(B), where HA is a Hilbert space associated with the first quantum system, and H_(B) is a Hilbert space associated with the second quantum system. The kets

represent the orthonormal eigenstates of the Hilbert space H_(A), and the kets

represents the orthonormal eigenstates of the Hilbert space H_(B), where i and j are positive integers. Any linear superposition of states in the Hilbert space H_(A){circle around (×)}H_(B) is given by: $\left. \Psi \right\rangle_{AB} = {\sum\limits_{i,j}{c_{ij}\left. i \right\rangle_{A}\left. j \right\rangle_{B}}}$ where the amplitudes c_(ij) are complex numbers satisfying the condition: ${\sum\limits_{ij}{c_{ij}}^{2}} = 1$

Special kinds of linear superpositions of states

are called “direct product states” and are represented by the product: $\left. \Psi \right\rangle_{AB} = {{\left. \psi \right\rangle_{A}\left. \psi \right\rangle_{B}} = {\left( {\sum\limits_{i}{c_{i}^{(A)}\left. i \right\rangle_{A}}} \right)\left( {\sum\limits_{j}{c_{j}^{(B)}\left. j \right\rangle_{B}}} \right)}}$ where

is a normalized linear superposition of states in the Hilbert space H_(A), and

is a normalized linear superposition of states in the Hilbert space H_(B). However, linear superpositions in the Hilbert space H_(A){circle around (×)}H_(B) that cannot be written as a product state are entangled states. In general, for a Hilbert space comprising two or more quantum subsystems, an entangled state is a linear superposition of states that cannot be written as a direct product state. The Bell states

and

are examples of entangled states, because the Bell states cannot be factored into products of the qubits

and

for any choice of the parameters α₁, β₁, α₂, and β₂.

Embodiments of the Present Invention

Various embodiments of the present invention are directed to generating entangled states of non-interacting qubits via interaction with a coherent state. FIG. 11 illustrates a control-flow diagram that describes one of many embodiments of the present invention for generating an entangled state of non-interacting qubits. The qubits can be photon-based qubits or matter-based qubits. For example, the basis states

and

can represent horizontal and vertical polarization states of an electromagnetic field, or the basis states

and

can represent the spin-up and spin-down states,

and

of an electron. In step 1102, an input state: $\left. \Psi \right\rangle_{in} = {\prod\limits_{n = 1}^{N}{\left. \psi \right\rangle_{n}\left. \alpha \right\rangle_{a}}}$ is generated, where $\left. \Psi \right\rangle_{n} = {\frac{1}{2^{1/2}}\left( {\left. 0 \right\rangle_{n} + \left. 1 \right\rangle_{n}} \right)}$ represents a qubit,

n is a qubit index, and

a identifies the field mode of the probe. The input state is a tensor product of N non-interacting qubits and a coherent state, which is called a “probe.” The input state can be rewritten as a linear superposition of states: $\left. \Psi \right\rangle_{in} = {\frac{1}{2^{N/2}}{\sum\limits_{j = 1}^{N + 1}{\left. \Gamma \right\rangle_{j}\left. \alpha \right\rangle_{a}}}}$ where $\left. \Gamma \right\rangle_{j} = {\sum\limits_{{({y_{1},\ldots\quad,y_{N}})} \in P_{j}}\left. {y_{1},\ldots\quad,y_{N}} \right\rangle_{j}}$ The components y_(n) of the N-tuple (y₁, . . . , y_(N)) are elements of the set {0,1} and corresponds to a tensor product |y₁, . . . ,

of the basis states

and

and P_(j) is a set of permutations associated with the N-tuple (y₁, . . . , y_(N)). The number of N-tuples in the set P_(j) is equal to the binomial coefficient $\begin{pmatrix} N \\ N_{1} \end{pmatrix},$ where N₁ is the number of “1” bits in the N-tuple (y₁, . . . y_(N)) For example, the set of permutations associated with the 3-tuple (1,0,0) is {(1,0,0),(0,1,0),(0,0,1)}. An example of an input state comprising 4 qubits is: $\left. \Psi \right\rangle_{in} = {\frac{1}{2^{4/2}}\left( {\left. 0 \right\rangle_{1} + \left. 1 \right\rangle_{1}} \right)\left( {\left. 0 \right\rangle_{2} + \left. 1 \right\rangle_{2}} \right)\left( {\left. 0 \right\rangle_{3} + \left. 1 \right\rangle_{3}} \right)\left( {\left. 0 \right\rangle_{4} + \left. 1 \right\rangle_{4}} \right)\left. \alpha \right\rangle_{a}}$ which can be rewritten as a linear superposition of states: $\left. \Psi \right\rangle_{in} = {\frac{1}{4}{\sum\limits_{j = 1}^{5}{\left. \Gamma \right\rangle_{j}\left. \alpha \right\rangle_{a}}}}$ where Γ⟩₁ = 0000⟩₁ Γ⟩₂ = 0000⟩₂ + 0100⟩₂ + 0010⟩₂ + 0001⟩₂ Γ⟩₃ = 1100⟩₃ + 1010⟩₃ + 1001⟩₃ + 0110⟩₃ + 0101⟩₃ + 0011⟩₃ Γ⟩₄ = 1110⟩₄ + 1101⟩₄ + 1011⟩₄ + 0111⟩₄ Γ⟩₅ = 1111⟩₅ In step 1104, the probe

is coupled to the product states |y₁, . . . ,

and the interaction is characterized by a quantum-mechanical, time-evolution operator: $\hat{U} = {\prod\limits_{n = 1}^{N}{\exp\left( \frac{{- {\mathbb{i}}}\quad{\hat{H}}_{a,n}t}{\hslash} \right)}}$ where

Ĥ_(a,n) is an interaction Hamiltonian that couples the probe

with the basis states

and

t is the interaction time. The interaction leaves the product states |y₁, . . . ,

unchanged. However, the probe experiences a phase shift or a translation represented by

as described below with reference to FIGS. 12A-12B. A phase shift is a change in the phase φ of the complex number α. A translation can be both a change in the average number of photons and a phase shift that is represented in phase space by adding a complex-valued parameter γ_(j) to the parameter α. The interaction between a probe and a product state is mathematically represented by: ${\hat{U}\left. {y_{1},\ldots\quad,y_{N}} \right\rangle_{j}\left. \alpha \right\rangle_{a}} = {{\sum\limits_{n = 1}^{N}\quad{{\exp\left( \frac{{- {\mathbb{i}}}{\hat{H}}_{n,a}t}{\hslash} \right)}\left. y_{n} \right\rangle_{j}\left. \alpha \right\rangle_{a}}} = {\left. {y_{1},\ldots\quad,y_{N}} \right\rangle_{j}\left. \alpha_{j} \right\rangle_{a}}}$ In step 1106, an output state

is obtained. The output state

is a linear superposition of states resulting from interaction of the probe with the product states of the input state

and is determined mathematically by applying the time-evolution Û operator to the input state

as follows: $\left. \Psi \right\rangle_{out} = {{\hat{U}\left. \Psi \right\rangle_{in}} = {{\frac{1}{2^{\frac{N}{2}}}{\sum\limits_{j = 1}^{N + 1}\quad{\hat{U}\left. \Gamma \right\rangle_{j}\left. \alpha \right\rangle_{a}}}} = {{\frac{1}{2^{\frac{N}{2}}}{\sum\limits_{j = 1}^{N + 1}\quad{\sum\limits_{{({y_{1},\ldots\quad,y_{N}})} \in P_{j}}\quad{\hat{U}\left. {y_{1},\ldots\quad,y_{N}} \right\rangle_{j}\left. \alpha \right\rangle_{a}}}}} = {{\frac{1}{2^{\frac{N}{2}}}{\sum\limits_{j = 1}^{N + 1}\quad{\sum\limits_{{({y_{1},\ldots\quad,y_{N}})} \in P_{j}}\quad{\left. {y_{1},\ldots\quad,y_{N}} \right\rangle_{j}\left. \alpha_{j} \right\rangle_{a}}}}} = {\frac{1}{2^{\frac{N}{2}}}{\sum\limits_{j = 1}^{N + 1}\quad{\left. \Gamma \right\rangle_{j}\left. \alpha_{j} \right\rangle_{a}}}}}}}}$ In step 1108, when the interaction results in a phase shift of the probe, the method proceeds to step 1110, otherwise when the interaction result in a translation of the probe, the method proceeds to step 1112. In step 1110, homodyne detection is performed on the probe

using an X quadrature measurement. The X quadrature measurement projects the output state

into an entangled state

which is a linear superposition of two or more states

In step 1114, the entangled state

is output. In step 1112, photon number detection is performed on the probe

The measurement projects the output state

into an entangled state

which is a linear superposition of two or more states

In step 1116, the entangled state

is output.

FIGS. 12A-12B are phase-space diagrams that represent a phase shift and a translation of a probe, respectively, as described above with reference to step 1104, in FIG. 11. In FIGS. 12A-12B, each probe is represented in phase space by a vector and an uncertainty disks, as described above with reference to FIG. 9. In FIG. 12A, a phase shift interaction between the probe

1202 and the product state |y₁, . . . ,

causes a phase shift g_(j)(θ) 1206 that results in a phase-shifted probe

1204. The average number of photons |α| associated with the probe

and with the phase-shifted probe

are identical, as represented by the identical length vectors 1202 and 1204. In FIG. 12B, a translation interaction between the probe

1202 and the product state |y₁, . . . ,

causes a translation γ_(j) 1206 that results in a translated probe |α+

1208. The translated probe

1208 has a phase angle φ′ and average number of photons |α+γ_(j)| that are different from the phase angle φ and the average number of photons |α| of the probe

1202.

The probe

and a reference coherent state

are initially prepared with identical average photon numbers and phase angles. After the probe interaction in step 1104, the phase shift or translation is determined by comparing the probe

to the reference coherent state

The probe

parameter a can be thought of as initially having a phase angle equal to zero. In other words, α is initially real valued and lies on the X quadrature axis described above with reference to FIG. 9.

FIG. 13 illustrates a first implementation 1300 of the method described above with reference to FIG. 11 for the 4-qubit input state

that represents an embodiment of the present invention. The implementation 1300 comprises an interaction region 1302 and homodyne detection 1304. The interaction region 1302 comprises four separate and identical interaction mediums 1307-1310 for interacting basis states of matter-based qubits with the probe. The interactions are characterized by an interaction Hamiltonian: Ĥ_(a,n) ¹=hχâ^(†){circumflex over (aσ)}_(z,n) where â^(†)â is the number operator of the probe

and χ is a constant that represents the coupling strength between the probe

and the basis states. The operator {circumflex over (σ)}_(z,n)=

−

is an inversion operator that operates on the basis states as follows: {circumflex over (σ)}_(z,n)

=

and {circumflex over (σ)}_(z,n)

The interaction Hamiltonian Ĥ_(a,n) ¹ arises from an electric or magnetic dipole interaction between the basis states of a matter-based qubit and the probe

The transmission channels 1311-1315 transmit the probe

and the qubits

and

into, and out of, the interaction region 1302. The transmission channels 1311-1315 are separate and prevent the basis states from interacting with one another. The probe

is transmitted to each of the interaction mediums 1307-1310 and interacts separately with each of the qubits, for a period of time t. The time-evolution operator that characterizes the interaction region 1302 is: ${\hat{U}}_{1} = {\prod\limits_{n = 1}^{4}\quad{\exp\left( {{- {\mathbb{i}\theta}}\quad{\hat{a}}^{\dagger}\hat{a}\quad{\hat{\sigma}}_{z,n}} \right)}}$ where θ is the interaction strength and is equal to the product χt. The state output from the interaction region 1302 is mathematically represented by applying the time-evolution operator Û₁ to the input state

The basis states

and

in each term of

are unchanged by interactions with the probe

However, the probe

accumulates a phase shift. For example, the product state

interaction with the probe

is: ${\prod\limits_{n = 1}^{4}\quad{{\exp\left( {{- {\mathbb{i}\theta}}\quad{\hat{a}}^{\dagger}\hat{a}\quad{\hat{\sigma}}_{z,n}} \right)}\left. 1000 \right\rangle\left. \alpha \right\rangle_{a}}} = {{{\exp\left( {{- {\mathbb{i}\theta}}\quad{\hat{a}}^{\dagger}\hat{a}\quad{\hat{\sigma}}_{z,1}} \right)}{\left. 1 \right\rangle_{1} \otimes {\exp\left( {{- {\mathbb{i}\theta}}\quad{\hat{a}}^{\dagger}\hat{a}\quad{\hat{\sigma}}_{z,2}} \right)}}{\left. 0 \right\rangle_{2} \otimes {\exp\left( {{- {\mathbb{i}\theta}}\quad{\hat{a}}^{\dagger}\hat{a}\quad{\hat{\sigma}}_{z,3}} \right)}}{\left. 0 \right\rangle_{3} \otimes {\exp\left( {{- {\mathbb{i}\theta}}\quad{\hat{a}}^{\dagger}\hat{a}\quad{\hat{\sigma}}_{z,4}} \right)}}\left. 0 \right\rangle_{4}\left. \alpha \right\rangle_{a}} = {\left. 1000 \right\rangle\left. {\alpha\quad{\mathbb{e}}^{- {\mathbb{i}2\theta}}} \right\rangle_{a}}}$ The interaction region 1302 outputs the linear superposition of states: $\left. \Psi \right\rangle_{out} = \quad{\frac{1}{4}\left\lbrack \quad{{\left. \Gamma \right\rangle_{1}\left. {\alpha\quad{\mathbb{e}}^{- {\mathbb{i}4\theta}}} \right\rangle_{a}} + {\left. \Gamma \right\rangle_{2}\left. {\alpha\quad{\mathbb{e}}^{- {\mathbb{i}2\theta}}} \right\rangle_{a}} + {\left. \Gamma \right\rangle_{3}\left. \alpha\quad \right\rangle_{a}} + {\left. \Gamma \right\rangle_{4}\left. {\alpha\quad{\mathbb{e}}^{- {\mathbb{i}2\theta}}} \right\rangle_{a}} + {\left. \Gamma \right\rangle_{5}\left. {\alpha\quad{\mathbb{e}}^{- {\mathbb{i}4\theta}}} \right\rangle_{a}}} \right\rbrack}$

FIG. 14A is a phase-space diagram of the probes in the output state

output from the interaction region 1302. In FIG. 14A, the vectors 1402-1406 represent the probes

respectively. In FIG. 13, the probe is output from the interaction region 1302 and input to the homodyne detection 1304. Operation of the homodyne detection 1304 is represented by projecting the vectors 1402-1406 and the associated uncertainty disks onto the X-quadrature axis 1408, as indicated by dashed-line directional arrows, such as directional arrow 1410. Uncertainty disks 1412-1414 identify ranges of X-quadrature values output from the homodyne detection 1304. An X-quadrature value that falls within one of the three uncertainty disks 1412-1414 reveals an associated entangled state of qubits output from the implementation 1400. For example, a homodyne detection 1304 that outputs an X-quadrature value within the range identified ay the uncertainty disk 1412 reveals the entangled state comprising a linear superposition of the states

and

A linear superposition of entangled states output from to the interaction region 1302 and input to the homodyne detection 1304 is: $\left. \Psi \right\rangle_{out}^{HD} = {{\frac{2^{\frac{1}{2}}{f\left( {x,{\alpha\quad\cos\quad 4\theta}} \right)}}{4}\left. \Sigma \right\rangle_{1}} + {\frac{{8^{\frac{1}{2}}{f\left( {x,{\alpha\quad\cos\quad 2\theta}} \right)}}\quad}{4}\quad\left. \Sigma \right\rangle_{2}} + {\frac{6^{\frac{1}{2}}{f\left( {x,\alpha} \right)}}{4}\left. \Sigma \right\rangle_{3}}}$ where ${{\left. \Sigma \right\rangle_{1} = {\frac{1}{2^{\frac{1}{2}}}\left( {{{\mathbb{e}}^{{\mathbb{i}\phi}_{1}{(ϰ)}}\left. 0000 \right\rangle} + {{\mathbb{e}}^{- {{\mathbb{i}\phi}_{1}{(ϰ)}}}\left. 1111 \right\rangle}} \right){\quad\quad}}}\quad{{\left. \Sigma \right\rangle_{2} = {\frac{1}{8^{\frac{1}{2}}}\left\lbrack {{{\mathbb{e}}^{{\mathbb{i}\phi}_{2}{(\theta)}}\left( {\left. 1000 \right\rangle + \left. 0100 \right\rangle + \left. 0010 \right\rangle + \left. 0001 \right\rangle} \right)} + {{\mathbb{e}}^{- {{\mathbb{i}\phi}_{2}{(\theta)}}}\left( {\left. 0111 \right\rangle + \left. 1011 \right\rangle + \left. 1101 \right\rangle + \left. 1110 \right\rangle} \right)}} \right\rbrack}},{and}}{\left. \Sigma \right\rangle_{3} = {\frac{1}{6^{\frac{1}{2}}}\left( {\left. 1100 \right\rangle + \left. 1010 \right\rangle + \left. 1001 \right\rangle + \left. 0110 \right\rangle + \left. 0101 \right\rangle + \left. 0011 \right\rangle} \right)}}}\quad$ are normalized entangled states output as a result of the homodyne detection 1304. The homodyne detection 1304 creates phase shifts in the product states of entangled states

and

where the phases are given by: φ₁(x)=α sin 4θ(x−2α cos 4θ) mod 2π, and φ₂(x)=α sin 2θ(x−2α cos 2θ) mod 2π Note that in FIG. 13, the phase shift created by the homodyne detection 1304 is represented by applying a phase shift e^(iφ(x)) to the qubit

where the phase φ(x) represents either the phase φ₁(x) or the phase φ₂ (x).

The phase shifts in the product states of the entangled states

and

can be corrected by applying single-qubit phase shift operations that are determined by the homodyne detection 1304 output. FIG. 15A illustrates a correction that removes the phase φ₁(x) from the entangled states output by the implementation 1300 are represents an embodiment of the present invention. In FIG. 15A, when the homodyne detection 1304 reveals the entangled state

the phase shift operation 1502 removes the phase φ₁(x). The phase shift operation 1502 is mathematically characterized by a phase shift operator: {circumflex over (Φ)}₁=exp(−iφ ₁(x){circumflex over (σ_(z,1)) )} The phase shift operation 1502 can be equally applied to any of the four qubits, because the final state is invariant under any permutation of the qubit labels. Application of the phase shift operation 1502 is mathematically represented by: ${{\Phi_{1}\left. \Sigma \right\rangle_{1}} = {\frac{1}{2^{\frac{1}{2}}}\left( {\left. 0000 \right\rangle + \left. 1111 \right\rangle} \right)\frac{8{{f\left( {ϰ,{\alpha\quad\cos\quad 2\theta}} \right)}}^{2}}{16}}},{{and}\quad\frac{2{{f\left( {ϰ,{\alpha\quad\cos\quad 4\theta}} \right)}}^{2}}{16}}$ FIG. 15B illustrates a correction that removes the phase φ₂(x) from the entangled states output by the implementation 1300 are represents an embodiment of the present invention. In FIG. 15B, when the homodyne detection 1306 reveals the entangled state

four phase shift operations 1504-1507 remove the phase φ₂(x). The phase shift operations 1504-1507 are mathematically characterized by a tensor product of phase shift operators: {circumflex over (Φ)}₁{circle around (×)}{circumflex over (Φ)}₂{circle around (×)}{circumflex over (Φ)}₃{circle around (×)}{circumflex over (Φ)}₄ where ${\hat{\Phi}}_{j} = {\exp\left( \frac{{- {{\mathbb{i}\phi}_{2}(ϰ)}}{\hat{\sigma}}_{z,j}}{2} \right)}$ The phase shift operations 1504-1507 break a single qubit operation into an operation that is applied to each of the four qubits. Application of the phase shift operations 1504-1507 is mathematically represented by: ${{{\hat{\Phi}}_{1} \otimes {\hat{\Phi}}_{2} \otimes {\hat{\Phi}}_{3} \otimes {\hat{\Phi}}_{4}}\left. \Sigma \right\rangle_{2}} = {\frac{1}{8^{\frac{1}{2}}}\left\lbrack \quad{\left( {\left. 1000 \right\rangle + \left. 0100 \right\rangle + \left. 0010 \right\rangle + \left. 0001 \right\rangle} \right) + \left( {\left. 0111 \right\rangle + \left. 1011 \right\rangle + \left. 1101 \right\rangle + \left. 1110 \right\rangle} \right)} \right\rbrack}$ Note that no phase shift correction is needed for the state: $\left. \sum \right\rangle_{3} = {\frac{1}{6^{1/2}}\left( {\left. 1100 \right\rangle + \left. 1010 \right\rangle + \left. 1001 \right\rangle + \left. 0110 \right\rangle + \left. 0101 \right\rangle + \left. 0011 \right\rangle} \right)}$ Alternatively, the phases in the final states can be noted and tracked through any system that uses the entangled state.

The amplitudes in the output state

are functions of the X-quadrature value x, and squaring the amplitudes gives probability distributions that when integrated over the X-quadrature axis reveal the probabilities associated with obtaining the entangled states

and

FIG. 14B is a plot of three probability distributions versus output values from the homodyne detection 1304 that represents one of many embodiments of the present invention. Curves 1416-1418 represent probability distributions: $\frac{6{{f\left( {x,\alpha} \right)}}^{2}}{16}$ respectively, that are obtained by squaring of the amplitudes of the entangled states |Σ>₁, |Σ>₂, and |Σ>₃ in |Ψ>_(out) ^(HD). The probability distributions 1416-1418 are centered about the X-quadrature values 2α, 2α cos 2θ, and 2α cos 4θ, respectively. The regions 1420-1422 correspond to the X-quadrature regions identified by the uncertainty disks 1414, 1413, and 1412, respectively, shown in FIG. 14A. Integrating each of the probability distributions 1416-1418 over the X-quadrature axis 1408 gives the corresponding probability that an X-quadrature value falls within the regions 1420-1422 and of obtaining the entangled states |Σ>₃, |Σ>₂, and |Σ>₁, respectively. The integral values of the probability distributions 1416-1418 indicate that there is ⅜ probability of outputting the entangled state |Σ>₃, a ½ probability of outputting the entangled state |Σ>₂, and a ⅛ probability of outputting the entangled state |Σ>₁, respectively.

FIG. 16 illustrates a second implementation 1600 of the method described above with reference to FIG. 11 for the 4-qubit input state |Ψ>_(in) that represents an embodiment of the present invention. The implementation 1600 comprises an interaction region 1602 and homodyne detection 1604. The interaction region 1602 comprises four separate and identical cross-Kerr interaction mediums 1607-1610 for interacting basis states of photon-based qubits with the probe and a rotate quadrature operation 1611. The cross-Kerr interaction mediums are characterized by an interaction Hamiltonian: Ĥ _(a,n) ² =hχâ ^(†) â({circumflex over (b)} ^(†) {circumflex over (b)})_(n) where ({circumflex over (b)}^(†){circumflex over (b)}) is a number operator that operates on the photon basis states as follows: ({circumflex over (b)} ^(†) {circumflex over (b)})_(n)|0>_(n)=0, and ({circumflex over (b)} ^(†) {circumflex over (b)})_(n)|1>_(n)=1·|1>_(n) The transmission channels 1612-1616 separately transmit the probe

and the photon qubits

and

into, and out of, the interaction region 1602. The transmission channels can be optical fibers that prevent the states from interacting. The probe

is transmitted to each of the Kerr interaction mediums 1607-1610 and interacts separately with each of the qubits, for a period of time t. The time-evolution operator characterizing the interactions is: ${\hat{U}}_{2} = {\prod\limits_{n = 1}^{4}\quad{\exp\left( {{- i}\quad\theta{\hat{\quad a}}^{\dagger}{\hat{a}\left( {{\hat{b}}^{\dagger}\hat{b}} \right)}_{n}} \right)}}$ Applying the time-evolution operator Û₂ to the input state

represents operation of the interaction mediums 1607-1610 and gives an output state comprising a linear superposition of states: $\begin{matrix} {\left. \Psi\quad \right\rangle_{out} = {\frac{1}{4}\left\lbrack {{\left. \Gamma\quad \right\rangle_{1}\left. \alpha \right\rangle_{a}} + {\left. \Gamma\quad \right\rangle_{2}\left. {\alpha\quad{\mathbb{e}}^{{- i}\quad\theta}} \right\rangle_{a}} +} \right.}} \\ \left. {{\left. \Gamma\quad \right\rangle_{3}\left. {\alpha\quad{\mathbb{e}}^{{- {i2}}\quad\theta}} \right\rangle_{a}} + {\left. \Gamma\quad \right\rangle_{4}\left. {\alpha\quad{\mathbb{e}}^{{- {i3}}\quad\theta}} \right\rangle_{a}} + {\left. \Gamma\quad \right\rangle_{5}\left. {\alpha\quad{\mathbb{e}}^{{- {i4}}\quad\theta}} \right\rangle_{a}}} \right\rbrack \end{matrix}$

FIG. 17A is a phase-space diagram of the probes in the output state

output from the interaction region 1602. In FIG. 17A, the vectors 1702-1706 represent to probes

and

respectively. In order to prepare the output state

for homodyne detection 1604, the rotate quadrature operation 1611 adjusts the phase of the probe. The rotate quadrature operation 1611 is characterized by the rotation operator exp (i2θâ^(†)â) and applying the rotation operator to the output state

out gives: $\begin{matrix} {\left. \Psi\quad \right\rangle_{out}^{rot} = {\frac{1}{4}\left\lbrack {{\left. \Gamma\quad \right\rangle_{1}\left. {\alpha\quad{\mathbb{e}}^{{i2}\quad\theta}} \right\rangle_{a}} + {\left. \Gamma\quad \right\rangle_{2}\left. {\alpha\quad{\mathbb{e}}^{i\quad\theta}} \right\rangle_{a}} +} \right.}} \\ \left. {{\left. \Gamma\quad \right\rangle_{3}\left. \alpha\quad \right\rangle_{a}} + {\left. \Gamma\quad \right\rangle_{4}\left. {\alpha\quad{\mathbb{e}}^{{- i}\quad\theta}} \right\rangle_{a}} + {\left. \Gamma\quad \right\rangle_{5}\left. {\alpha\quad{\mathbb{e}}^{{- i}\quad 2\theta}} \right\rangle_{a}}} \right\rbrack \end{matrix}$

FIG. 17B is a phase-space diagram of the output state

In FIG. 17B, the vectors 1708-1712 represent rotated probes

and

respectively. The homodyne detection 1604 projects the uncertainty disks associated with the probes onto the X-quadrature axis 1714, as indicated by dashed-line directional arrows. An X-quadrature value output from the homodyne detection 1604 falls within one of the ranges identified by three uncertainty disks 1715-1717 and reveals the entangled state of qubits output from the implementation 1600. For example, a homodyne detection 1604 that outputs an X-quadrature value within the range of values identified by the uncertainty disk 1716 results in an entangled state comprising a linear superposition of the states

and

A linear superposition of entangled states input to the homodyne detection 1604 is: $\begin{matrix} {\left. \Psi\quad \right\rangle_{out}^{{H\quad D}\quad} = {{\frac{2^{1/2}{f\left( {x,{\alpha\quad{\cos 2}\quad\theta}} \right)}}{4}\left. \sum \right\rangle_{1}} + {\frac{8^{1/2}{f\left( {x,{\alpha\quad\cos\quad\theta}} \right)}}{4}\left. \sum \right\rangle_{2}} +}} \\ {\frac{6^{1/2}{f\left( {x,\alpha} \right)}}{4}\left. \sum \right\rangle_{3}} \end{matrix}$ where $\left. \sum \right\rangle_{1} = {\frac{1}{2^{1/2}}\left( {{{\mathbb{e}}^{i\quad{\phi_{1}{(x)}}}\left. 0000 \right\rangle} + {{\mathbb{e}}^{{- i}\quad{\phi_{1}{(x)}}}\left. 1111 \right\rangle}} \right)}$ $\begin{matrix} {\left. \sum \right\rangle_{2} = {\frac{1}{8^{1/2}}\left\lbrack {{{\mathbb{e}}^{i\quad{\phi_{2}{(\theta)}}}\left( {\left. 1000 \right\rangle + \left. 0100 \right\rangle + \left. 0010 \right\rangle + \left. 0001 \right\rangle} \right)} +} \right.}} \\ {{{\mathbb{e}}^{{- i}\quad{\phi_{2}{(\theta)}}}\left( {\left. 0111 \right\rangle + \left. 1011 \right\rangle + \left. 1101 \right\rangle + \left. 1110 \right\rangle} \right)},{and}} \end{matrix}$ $\left. \sum \right\rangle_{3} = {\frac{1}{6^{1/2}}\left( {\left. 1100 \right\rangle + \left. 1010 \right\rangle + \left. 1001 \right\rangle + \left. 0110 \right\rangle + \left. 0101 \right\rangle + \left. 0011 \right\rangle} \right)}$ are normalized entangled states. The homodyne detection 1604 creates phase shifts in the product states of

and

where the phases are given by: φ₁(x)=α sin 2θ(x−2α cos 2θ) mod 2π, and φ₂(x)=α sin θ(x−2α cos θ) mod 2π

The phase shifts in the product states of the entangled states

and

can be corrected by applying single-qubit phase shift operations that are determined by the homodyne detection 1604 output. FIG. 18A illustrates a correction that removes the phase φ₁(x) from the entangled state

output by the implementation 1600 and represents an embodiment of the present invention. In FIG. 16A, when the homodyne detection 1604 reveals the entangled state

the phase shift operation 1802 removes the phase φ₁(x). The phase shift operation 1802 is mathematically characterized by a phase shift operator: {circumflex over (Φ)}₁=exp(−iφ ₁(x){circumflex over (σ)}_(z,1)) The phase shift operation 1802 can be equally applied to any of the four qubits, because the final state is invariant under any permutation of the qubit labels. Application of the phase shift operation 1802 is mathematically represented by: ${{\hat{\Phi}}_{1}\left. \sum \right\rangle_{1}} = {\frac{1}{2^{1/2}}\left( {\left. 0000 \right\rangle + \left. 1111 \right\rangle} \right)}$ FIG. 18B illustrates a correction that removes the phase φ₂(x) from the entangled states output by the implementation 1600 and represents an embodiment of the present invention. In FIG. 18B, when the homodyne detection 1604 reveals the entangled state

four phase shift operations 1804-1807 remove the phase φ₂ (x). The phase shift operations 1804-1807 are mathematically characterized by a tensor product of phase shift operators: {circumflex over (Φ)}₁{circle around (×)}{circumflex over (Φ)}₂{circle around (×)}{circumflex over (Φ)}₃{circle around (×)}{circumflex over (Φ)}₄ where ${\hat{\Phi}}_{j} = {\exp\left( \frac{{- i}\quad{\phi_{2}(x)}{\hat{\sigma}}_{z,j}}{2} \right)}$ The phase shift operations 1804-1807 break a single qubit operation into an operation that is applied to each of the four qubits. Application of the phase shift operations 1804-1807 is mathematically represented by: $\begin{matrix} {{{{\hat{\Phi}}_{1} \otimes {\hat{\Phi}}_{2} \otimes {\hat{\Phi}}_{3} \otimes {\hat{\Phi}}_{4}}\left. \sum \right\rangle_{2}} = {\frac{1}{8^{1/2}}\left\lbrack \left( {\left. 1000 \right\rangle + \left. 0100 \right\rangle + \left. 0010 \right\rangle +} \right. \right.}} \\ {\left. \left. 0001 \right\rangle \right) + \left( {\left. 0111 \right\rangle +} \right.} \\ \left. \left. {\left. 1011 \right\rangle + \left. 1101 \right\rangle + \left. 1110 \right\rangle} \right) \right\rbrack \end{matrix}$ Note that no phase shift correction is needed for the entangled state: $\left. \sum \right\rangle_{3} = {\frac{1}{6^{1/2}}\left( {\left. 1100 \right\rangle + \left. 1010 \right\rangle + \left. 1001 \right\rangle + \left. 0110 \right\rangle + \left. 0101 \right\rangle + \left. 0011 \right\rangle} \right)}$ Alternatively, the phases in the final states can be noted and tracked through any system that uses the entangled state.

FIG. 17C is a plot of three probability distributions versus output values from the homodyne detection 1604 that represents one of many embodiments of the present invention. Curves 1718-1720 are probability distributions that are obtained by squaring the amplitudes of the entangled states

and

in

and are centered about the X-quadrature values 2α, 2α cos 2θ, and 2α cos 4θ, respectively. The regions 1722-1724 correspond to the X-quadrature regions identified by the uncertainty disks 1715-1717, respectively, shown in FIG. 17C. Integrating each of the probability distributions 1718-1720 over the X-quadrature axis 1714 gives the corresponding probability of an X-quadrature value falling within one of the regions 1722-1724 and of obtaining one of the entangled states

and

respectively. The integral values of the probability distributions 1718-1720 indicates that there is ⅜ probability of outputting the entangled state

a ½ probability of outputting the entangled state

and a ⅛ probability of outputting the entangled state

FIG. 19 illustrates a third implementation 1900 of the method described above with reference to FIG. 11 for the 4-qubit input state

that represents an embodiment of the present invention. The implementation 1900 comprises an interaction region 1902 and a photon number detector 1904. The interaction region 1902 comprises four separate and identical interaction mediums 1906-1909 and a displacement operation 1910. The interaction mediums are characterized by an interaction Hamiltonian: Ĥ _(a,n) ^(φ) =hχ{circumflex over (σ)} _(z,n)(â ^(†) e ^(iφ) +âe ^(−iφ)) where â^(†)e^(iφ)+âe^(−iφ) is a quadrature operator of the probe and is determined by the angle φ. When φ equals −π/2, the quadrature operator is the Y-quadrature operator, and when φ equals 0, the quadrature operator is the X-quadrature operator. The interaction Hamiltonian Ĥ_(a,n) ^(φ) represents a dipole coupling of a matter-based qubit with the probe. The transmission channels 1912-1916 transmit the probe

and the qubits

and

into, and out of, the interaction region 1902. The probe

is transmitted to each of the interaction mediums 1906-1909 and interacts separately with each of the qubits, for a period of time t. Interactions in the interaction region 1902 are characterized by the time-evolution operator: $\begin{matrix} {{\hat{U}}_{\varphi} = {\prod\limits_{n = 1}^{4}\quad{\exp\left( \frac{{- {\mathbb{i}}}{\hat{H}}_{a,n}^{\varphi}t}{\hslash} \right)}}} \\ {= {\prod\limits_{n = 1}^{4}{\hat{D}\left( {{\hat{\sigma}}_{z,n}\gamma\quad(\varphi)} \right)}}} \end{matrix}$ where {circumflex over (D)}({circumflex over (σ)}_(z,n)γ(φ))=exp({circumflex over (σ)}_(z,n)γ(φ)â ^(†)−{circumflex over (σ)}_(z,n)γ*(φ)â) is called the “displacement operator,” and γ(φ)=θe ^(i(φ−*/2)).

In an embodiment of the present invention, the angle φ is equal to −π/2. Applying the time-evolution operator Û_(−*/2) to the input state

represents the operation performed by the interaction region 1902 and gives an output state comprising a linear superposition of states: $\begin{matrix} {\left. \Psi \right\rangle_{out} = {{\hat{U}}_{{- \pi}/2}\left. \Psi \right\rangle_{in}}} \\ {= {\frac{1}{4}\left\lbrack {{\left. \Gamma \right\rangle_{1}\left. {\alpha + {4\theta}} \right\rangle_{a}} + {\left. \Gamma \right\rangle_{2}\left. {\alpha + {2\theta}} \right\rangle_{a}} + {\left. \Gamma \right\rangle_{3}\left. \alpha \right\rangle_{a}} +} \right.}} \\ \left. {{\left. \Gamma \right\rangle_{4}\left. {\alpha - {2\theta}} \right\rangle_{a}} + {\left. \Gamma \right\rangle_{5}\left. {\alpha - {4\theta}} \right\rangle_{a}}} \right\rbrack \end{matrix}$ FIG. 20A is a phase-space diagram of the output state

In FIG. 20A, the uncertainty disks 2002-2006 located along the X-quadrature axis correspond to the probes

and

respectively. The probes are each located at a different distance from the origin, which indicates that each of the probes corresponds to a different number of photons.

The states

and

in the output state

are not entangled states. In order to project onto a linear superposition of the states

and

the value of a in the probes states of the output state |Ψ>_(out) is set to zero. In FIG. 19, the displacement operation 1910 is applied to the probe in order to remove the average number of photons associated with the parameter α. The displacement operation 1910 is mathematically represented by the displacement operator: {circumflex over (D)}(−α)={circumflex over (D)}*(α)=exp(−αâ+α*â) Applying the displacement operator {circumflex over (D)}(−α) to the output state

gives a phase adjusted output state as follows: $\begin{matrix} {\left. \Psi \right\rangle_{out}^{adj} = {{\hat{D}\left( {- \alpha} \right)}\left. \Psi \right\rangle_{out}}} \\ {= {\frac{1}{4}\left\lbrack {{\left. \Gamma \right\rangle_{1}{\hat{D}\left( {- \alpha} \right)}\left. {\alpha + {4\theta}} \right\rangle_{a}} + {\left. \Gamma \right\rangle_{2}{\hat{D}\left( {- \alpha} \right)}\left. {\alpha + {2\theta}} \right\rangle_{a}} +} \right.}} \\ {{\left. \Gamma \right\rangle_{3}{\hat{D}\left( {- \alpha} \right)}\left. \alpha \right\rangle_{a}} + {\left. \Gamma \right\rangle_{4}{\hat{D}\left( {- \alpha} \right)}\left. {\alpha - {2\theta}} \right\rangle_{a}} +} \\ \left. {\left. \Gamma \right\rangle_{5}{\hat{D}\left( {- \alpha} \right)}\left. {\alpha - {4\theta}} \right\rangle_{a}} \right\rbrack \\ {= {\frac{1}{4}\left\lbrack {{\left. \Gamma \right\rangle_{1}\left. {4\theta} \right\rangle_{a}} + {\left. \Gamma \right\rangle_{2}\left. {2\theta} \right\rangle_{a}} + {\left. \Gamma \right\rangle_{3}\left. 0 \right\rangle_{a}} + {\left. \Gamma \right\rangle_{4}\left. {{- 2}\theta} \right\rangle_{a}} +} \right.}} \\ \left. {\left. \Gamma \right\rangle_{5}\left. {{- 4}\theta} \right\rangle_{a}} \right\rbrack \end{matrix}$

FIG. 20B is a phase-space diagram of the phase-adjusted output state

As a result of the displacement operation 1910, the probe

is translated to the vacuum state

and is located at the origin 2008. The probes

and |α−

in FIG. 20A, are translated to

and

respectively, which are at identical distances from the origin and have identical numbers of photons as indicated by circle 2010. The probes

and

in FIG. 20A, are translated to

and

which are also at identical distances from the origin and have identical numbers of photons as indicated by circle 2012.

In FIG. 19, photon number detector 1904 is used to determine the number of photons in the probe, and, as a result, the phase-adjusted output state

is projected into one of the number states represented in FIG. 20B. The photon number detector 1904 projects the phase adjusted output state

into: $\left. \Psi \right\rangle_{out}^{PNC} = {{\frac{1}{2^{3/2}}\left. \Phi \right\rangle_{1}} + {\frac{1}{2^{1/2}}\left. \Phi \right\rangle_{2}} + {\left( \frac{3}{2} \right)^{1/2}\left. \Phi \right\rangle_{3}}}$ where $\left. \Phi \right\rangle_{1} = {\frac{1}{2^{1/2}}\left( {\left. 0000 \right\rangle + {\left( {- 1} \right)^{m}\left. 1111 \right\rangle}} \right)}$ $\begin{matrix} {\left. \Phi \right\rangle_{2} = {\frac{1}{8^{1/2}}\left\lbrack {\left. 1000 \right\rangle + \left. 0100 \right\rangle + \left. 0010 \right\rangle + \left. 0001 \right\rangle + {\left( {- 1} \right)^{m}\left( {\left. 1110 \right\rangle +} \right.}} \right.}} \\ \left. \left. {\left. 1101 \right\rangle + \left. 1011 \right\rangle + \left. 0111 \right\rangle} \right) \right\rbrack \end{matrix}$ $\quad{\left. \Phi \right\rangle_{3} = {\frac{1}{6^{1/2}}\left( {\left. 1100 \right\rangle + \left. 1010 \right\rangle + \left. 1001 \right\rangle + \left. 0110 \right\rangle + \left. 0101 \right\rangle + \left. 0011 \right\rangle} \right)}}$ are the normalized entangled states, and m is the number of photons measured.

The phase shift (−1)^(m) in the entangled states

and

can be corrected by applying single-qubit phase shift operations. FIG. 21A illustrates a correction that removes the phase (−1)^(m) from the entangled state

output by the implementation 1900 and represents an embodiment of the present invention. In FIG. 21A, when the photon number detector 1904 reveals the state

the phase shift operation 2102 removes the phase (−1)^(m). The phase shift operation 2102 is mathematically characterized by a phase shift operator: ${\hat{\Omega}}_{1} = {{\exp\left( \frac{{- {\mathbb{i}m}}\quad\pi}{2} \right)}{\exp\left( \frac{{\mathbb{i}m}\quad\pi\quad{\hat{\sigma}}_{z,1}}{2} \right)}}$ The phase shift operation 2102 can be equally applied to any of the four qubits, because the final state is invariant under any permutation of the qubit labels. Application of the phase shift operation 2102 is mathematically represented by: ${{\hat{\Omega}}_{1}\left. \Phi \right\rangle_{1}} = {\frac{1}{2^{1/2}}\left( {\left. 0000 \right\rangle + \left. 1111 \right\rangle} \right)}$ FIG. 21B illustrates a correction that removes the phase (−1)^(m) from the entangled state

output by the implementation 1900 and represents an embodiment of the present invention. In FIG. 21B, when the photon number detector 1904 reveals the entangled state

the phase shift operations 2104-2107 remove the phase φ₂ (x). The phase shift operations 2104-2107 are mathematically characterized by a tensor product of phase shift operators: ${\exp\left( \frac{{- {\mathbb{i}m}}\quad\pi}{2} \right)}{{\hat{\Omega}}_{1} \otimes {\hat{\Omega}}_{2} \otimes {\hat{\Omega}}_{3} \otimes {\hat{\Omega}}_{4}}$ where ${\hat{\Omega}}_{j} = {\exp\left( \frac{{\mathbb{i}m}\quad\pi\quad{\hat{\sigma}}_{z,j}}{4} \right)}$ The phase shift operations 2104-2107 break a single qubit operation into an operation that is applied to each of the four qubits. The phase shift operations 2104-2107 applied to the entangled state

is mathematically represented by: $\left. {{{\exp\left( \frac{{- {\mathbb{i}}}\quad m\quad\pi}{2} \right)}{{\hat{\Omega}}_{1} \otimes {\hat{\Omega}}_{2} \otimes {\hat{\Omega}}_{3} \otimes {\hat{\Omega}}_{4}}}❘\Phi} \right\rangle_{2} = {\frac{1}{8^{1/2}}\begin{bmatrix} {\left. \left. {\left. {{{\left. \left( {1000} \right. \right\rangle +}}0100} \right\rangle + \left. 0010 \right\rangle + {{\quad }0001}} \right\rangle \right) +} \\ \left. \left. {{{\left. {{{\left. {{{\left( {0111} \right\rangle +}}1011} \right\rangle +}}1101} \right\rangle +}}1110} \right\rangle \right) \end{bmatrix}}$ Note that no phase shift correction is needed for the state: $\left. \left. {{{\left. {{{\left. {{{\left. {{{\left. {\left. {❘\Sigma}\quad \right\rangle_{3} = {\frac{1}{6^{1/2}}\left( {\left. 1100 \right\rangle +} \right.1010}} \right\rangle +}}1001} \right\rangle +}}0110} \right\rangle +}}0101} \right\rangle +}}0011} \right\rangle \right)$ Alternatively, the phases in the final states can be noted and tracked through any system that uses the entangled state. Note that the photon number detector 1904 absorbs the photons, and outputs a signal that can be processed by digital circuits. Photomultiplier tubes, avalanche photodiodes, and high efficiency photon detection devices, such as a visible light photon counter, can be used to determine the number of photons in the probe number states of the output state

Squaring the amplitudes of the output state

gives the probability of outputting the entangled states

and

For example, there is a ⅛ probability of outputting the entangled state

a ½ probability of outputting the entangled state

and a ⅜ probability of outputting the entangled state

In an embodiment of the present invention, the angle φ is equal to zero. Applying the time-evolution operator Û₀ to the input state

represents operation of the interaction region 1902 and gives an output state comprising the linear superposition of states: $\begin{matrix} \left. {\left. {❘\Psi} \right\rangle_{out} = {{\hat{U}}_{0}❘\Psi}} \right\rangle_{in} \\ {= {\frac{1}{4}\begin{bmatrix} {\left. {{{\left. {{{\left. {{{\left. {{{\left. \Gamma \right\rangle_{1}}\alpha} + {{\mathbb{i}}\quad 4\theta}}\quad \right\rangle_{a} +}}\Gamma} \right\rangle_{2}}\alpha} + {{\mathbb{i}}\quad 2\theta}} \right\rangle_{a} + \left. \Gamma \right\rangle_{3}}}\alpha} \right\rangle_{a} +} \\ \left. \quad{{{\quad{{\left. \Gamma \right\rangle_{4}\left. {\alpha - {{\mathbb{i}}\quad 2\theta}} \right\rangle_{a}} + \left. \Gamma \right\rangle_{5}}}\alpha} - {{\mathbb{i}}\quad 4\theta}} \right\rangle_{a} \end{bmatrix}}} \end{matrix}$ FIG. 22A is a phase-space diagram of the output state

In FIG. 22A, the uncertainty disks 2202-206 are aligned vertically above and below the X-quadrature axis and correspond to the probes

and

respectively.

The states

and

are entangled by projecting the output state

onto a linear superposition of the states

and

which is achieved by translating the probes to the Y-quadrature axis. In FIG. 19, the displacement operation 1910 translates the probes in

to the Y-quadrature axis as follows: $\begin{matrix} \left. {\left. {❘\Psi} \right\rangle_{out}^{adj} = {{\hat{D}\left( {- \alpha} \right)}❘\Psi}} \right\rangle_{in} \\ {= {\frac{1}{4}\begin{bmatrix} {\left. {{{\left. {{{\left. {{\left. \Gamma \right\rangle_{1}}{\mathbb{i}}\quad 4\theta} \right\rangle_{a} + \left. \Gamma \right\rangle_{2}}}{\mathbb{i}}\quad 2\theta} \right\rangle_{a} + \left. \Gamma \right\rangle_{3}}}0} \right\rangle_{a} +} \\ \left. {{{\left. {{\left. \Gamma \right\rangle_{4}} - {{\mathbb{i}}\quad 2\theta}} \right\rangle_{a} + \left. \Gamma \right\rangle_{5}}} - {{\mathbb{i}}\quad 4\theta}} \right\rangle_{a} \end{bmatrix}}} \end{matrix}$

FIG. 22B is a phase-space diagram of the phase-adjusted output state

Circle 2208 represents the probes

and

which are both the same distance form the origin and have identical numbers of photons. Circle 2210 represents the probes

and

which are both the same distance form the origin and have identical numbers of photons. The vacuum state

is located at the origin 2212.

In FIG. 19, photon number detector 1904 gives the number of photons in the probe, which projects the phase-adjusted output state

into one of the probes represented in FIG. 22B. The state output by the photon number detector 1904 is: $\left. {{{\left. {{{\left. {❘\Psi} \right\rangle_{out}^{PNC} = \frac{1}{2^{3/2}}}}\Phi} \right\rangle_{1} + {\frac{1}{2^{1/2}}\left. \Phi \right\rangle_{2}} + \left( \frac{3}{2} \right)^{1/2}}}\Phi} \right\rangle_{3}$ which is identical to the output state described above for phase angle φ equal to −π/2.

Although the present invention has been described in terms of particular embodiments, it is not intended that the invention be limited to these embodiments. Modifications within the spirit of the invention will be apparent to those skilled in the art. For example, in an alternate embodiment of the present invention, rather than employing interaction regions having two or more interaction mediums as described above with reference to FIGS. 13, 16, and 19, the interaction region can have a single interaction medium. FIG. 23 illustrates an interaction region comprising a single interaction medium that represents one of many embodiments of the present invention. In FIG. 23, the probe

is repeatedly transmitted into an interaction medium 2302 via a transmission channel 2304. Transmission of the probe into the interaction medium 2302 is synchronized with the transmission of the qubits

and

so that the probe interacts separately with each qubit. In an alternate embodiment of the present invention, rather than performing the displacement operation 1910, in FIG. 19, on the probe after the probe is output from the interaction region 1902, the displacement operation can be performed before the probe is input to the interaction region 1902. In an alternate embodiment of the present invention, the displacement operation 1904 can be eliminated by using the vacuum state with α equal to “0.” In an alternate embodiment of the present invention, the use of other interaction mediums characterized by different interaction Hamiltonians can be used. In alternate embodiments of the present invention, multiple probe beams can be used to interact with subsets of qubits and then interfere at a beam-splitter prior to a measurement. In alternate embodiments of the present invention, the entangled qubits can be quantum dots, magnetic moments, such as superconducting flux or magnetic spin qubits, or the qubits can be modes of nanoelectromechanical systems. In alternate embodiments of the present invention, each of the N matter-based qubits can be trapped within one of the N separate interaction mediums in the interaction region. The probe is transmitted to each of the N interaction mediums to separately interact with the N qubits. In an alternate embodiment of the present invention, an interaction region comprises numerous different kinds of interaction mediums. For example, an interaction region comprises a combination of interaction mediums for coupling matter-based qubits to a probe and Kerr interaction mediums for coupling photon-based qubits to the probe.

The foregoing description, for purposes of explanation, used specific nomenclature to provide a thorough understanding of the invention. However, it will be apparent to one skilled in the art that the specific details are not required in order to practice the invention. The foregoing descriptions of specific embodiments of the present invention are presented for purposes of illustration and description. They are not intended to be exhaustive of or to limit the invention to the precise forms disclosed. Obviously, many modifications and variations are possible in view of the above teachings. The embodiments are shown and described in order to best explain the principles of the invention and its practical applications, to thereby enable others skilled in the art to best utilize the invention and various embodiments with various modifications as are suited to the particular use contemplated. It is intended that the scope of the invention be defined by the following claims and their equivalents: 

1. A method for preparing an entangled state of qubits, the method comprising: providing a probe and N non-interacting qubits, each qubit comprises a linear superposition of two basis states; transmitting the probe into an interaction region that separately couples the probe to each of the qubits and produces a number of different probe states; outputting from the interaction region a linear superposition of states, each state in the linear superposition of states comprises a tensor product of entangled basis states and one of the different probe states; and projecting the linear superposition of states into one of the entangled states by measuring the state of the probe.
 2. The method of claim 1 wherein the probe is a coherent state and the N non-interacting qubits further comprises one or more of: matter-based qubits; photon-based qubits; a linear superpositions of two quantum dot states; and a linear superposition of two magnetic states.
 3. The method of claim 1 wherein the detector is a photon number detector.
 4. The method of claim 1 wherein the detector is an X quadrature homodyne detector.
 5. The method of claim 1 wherein the interaction region further comprises one or more dipole coupling mediums, each medium characterized by an interaction Hamiltonian: Ĥ=hχâ^(†){circumflex over (aσ)}_(z) where χ is a constant that represents the coupling strength between the probe and basis states of the qubit, â^(†)â is a number operator that operates on the probe, and


6. The method of claim 1 wherein the interaction region further comprises one or more cross-Kerr interaction mediums, each medium characterized by a cross-Kerr interaction Hamiltonian: Ĥ=hχâ^(†)âb^(†){circumflex over (b)} where χ is a constant that represents the coupling strength between the probe and basis states of the qubit, â^(†)â is a number operator that operates on the probe, and {circumflex over (b)}^(†){circumflex over (b)} is a number operator that operates on a photon basis state.
 7. The method of claim 1 wherein the interaction region further comprises one or more interaction mediums, each medium characterized by an interaction Hamiltonian: Ĥ=hχ{circumflex over (σ)} _(z)(â ^(†) e ^(iφ) +âe ^(−iφ)) where χ is a constant that represents the coupling strength between the probe and basis states of the qubit,

and â^(†)e^(iφ)+âe^(−iφ) is a quadrature operator of the probe and is determined by the angle φ.
 8. The method of claim 1, wherein the method further comprises a correction operation that removes phase shifts acquired by the entangled states.
 9. A method for preparing an entangled state of qubits, the method comprising: providing a probe and N non-interacting qubits, each qubit comprises a linear superposition of two basis states; transmitting the probe into an interaction region that separately couples the probe to each of the qubits and produces a number of different phase-shifted probe states; outputting from the interaction region a linear superposition of states, each state in the linear superposition of states comprises a tensor product of a linear superposition of entangled states and one of the phase-shifted probe states; and transmitting the phase-shifted probe state to a homodyne detector that projects the linear superposition of states into one of the entangled states by measuring the state of the coherent state.
 10. The method of claim 9, wherein the probe is a coherent state and the N non-interacting qubits further comprises one or more of: matter-based qubits; photon-based qubits; a linear superpositions of two quantum dot states; and a linear superposition of spin magnetic moments.
 11. The method of claim 9, wherein the detector is an X quadrature homodyne detector.
 12. The method of claim 9, wherein the interaction region further comprises dipole coupling mediums, each medium characterized by an interaction Hamiltonian: Ĥ=hχâ^(†){circumflex over (aσ)}_(z) where â^(†)â is a number operator that operates on the probe, χ is a constant that represents the coupling strength between the probe and basis states of the qubit, and


13. The method of claim 9, wherein the interaction region further comprises cross-Kerr interaction mediums, each medium characterized by a cross-Kerr Hamiltonian: Ĥ=hχâ^(†){circumflex over (ab)}^(†){circumflex over (b)} where χ is a constant that represents the coupling strength between the probe and basis states of the qubit, â^(†)â is a number operator that operates on the probe, and {circumflex over (b)}^(†){circumflex over (b)} is a number operator that operates on a photon basis state.
 14. The method of claim 13, wherein the interaction region further comprises applying a quadrature rotation operation to the probe output from the interaction region.
 15. The method of claim 9, wherein the method further comprises a correction operation that removes phase shifts acquired by the entangled states.
 16. A method for preparing an entangled state of qubits, the method comprising: providing a probe and N non-interacting qubits, each qubit comprises a linear superposition of two basis states; transmitting the probe into an interaction region that separately couples the probe to each of the qubits and produces a number of different translated coherent states; outputting from the interaction region a linear superposition of states, each state in the linear superposition of states comprises a tensor product of a linear superposition of entangled states and one of the translated probe states; and transmitting the displaced probe state to a photon number detector that projects the linear superposition of states into one of the entangled states by measuring the state of the coherent state.
 17. The method of claim 16 wherein the probe is a coherent state and the N non-interacting qubits further comprise one or more of: matter-based qubits; a linear superpositions of two quantum dot states; and a linear superposition of spin magnetic moments.
 18. The method of claim 16 wherein the detector is a photon number detector.
 19. The method of claim 16 wherein the interaction region further comprises one or more interaction mediums, each medium characterized by an interaction Hamiltonian: Ĥ=hχ{circumflex over (σ)} _(z)(â ^(†) e ^(iφ) +âe ^(−iφ)) where χ is a constant that represents the coupling strength between the probe and basis states of the qubit,

and â^(†)e^(iφ)+âe^(−iφ) is a quadrature operator of the probe and is determined by the angle φ.
 20. The method of claim 19, wherein the interaction region further comprises applying a probe displacement operation to the probe output from the interaction region.
 21. The method of claim 16, wherein the method further comprises a correction operation that removes phase shifts acquired by the entangled states. 